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I .  INTRODUCTION  AND  BACKGROUND 

One  of  the  most  significant  advances  in  objective  weather 
prediction,  since  the  introduction  of  numerical  weather 
prediction  in  the  1950's  and  satellite  remote  sensing  capa¬ 
bilities  in  the  1960's,  has  been  the  development  of  Model 
Output  Statistics  (MOS)  weather  forecasting  method  by  Glahn 
and  Lowry  (1972).  In  general,  this  technique  is  the  deter¬ 
mination  of  a  statistical  relationship  between  an  operational 
weather  element  (predictand) ,  which  may  or  may  not  be  fore¬ 
cast  by  numerical  methods,  and  numerical  model  output  varia¬ 
bles  (predictors),  usually  via  linear  regression  methods. 

The  resulting  predictand/predictor  regression  equations 
provide  the  basis  for  generating  a  statistical  weather  pre¬ 
diction.  The  National  Weather  Service  (NWS)  has  included  MOS 
as  an  integral  part  of  their  weather  forecasting  operations 
since  the  early  1970's.  Currently,  the  NWS  maintains  MOS 
prediction  equations  for  approximately  15  weather  elements 
(e.g.  ceiling,  visibility,  obstructions  to  vision,  precipi¬ 
tation,  etc.)  at  forecast  times  ranging  from  6  to  48  hours. 
These  forecasts  are  routinely  provided  to  approximately  295 
civilian  and  190  military  locations  throughout  the  continental 
United  States  (CONUS)  and  Alaska  [Glahn,  1987] . 

Based  on  the  impressive  results  achieved  with  the  NWS 
MOS  program,  the  Department  of  Defense  (DOD) ,  through  the  Air 


Weather  Service  (AWS) ,  implemented  and  operated  a  quasi- 
global  version  of  the  NWS  MOS  system  at  the  Air  Force  Global 
Weather  Center  (AFGWC) ,  Offutt  AFB ,  Nebraska  [Best  and  Pryor, 
1983]  .  The  first  operational  forecasts  obtained  from  the  AWS 
MOS  system  were  produced  by  AFGWC  in  December  1980  and  the 
system  ran  operationally  for  a  period  of  approximately  18 
months.  Regions  for  which  operational  MOS  forecasts  were 
produced  included  Europe,  Asis  (including  Korea  and  Japan), 
the  South  China  Sea  (including  the  Philippines  and  Taiwan) , 
the  near  and  middle  east  and  northern  Africa.  The  AWS  MOS 
program  was  terminated  with  the  recent  decision  to  replace  the 
current  hemispheric  primitive  equation  (PE)  model  with  a 
spectral  global  dynamic  model  [Klein,  1981] . 

Throughout  its  tenure  as  an  operational  forecast  scheme, 
the  AWS  MOS  system  provided  the  U.S.  Air  Force  with  a  rela¬ 
tively  low  cost,  flexible  and  responsive  prediction  network. 
Further  development  of  the  AWS  MOS  system  has  been  postponed 
until  sufficient  spectral  model  output  is  archived. 

The  U.S.  Navy,  by  virtue  of  its  unique  marine  forecast¬ 
ing  responsibilities,  has  a  keen  interest  in  applying  MOS 
forecasting  schemes  to  global  oceanic  regions.  Through  the 
research  and  development  efforts  of  the  Naval  Environmental 
Prediction  Research  Facility  (NEPRF)  in  Monterey,  California, 
the  Navy  has  sponsored  a  limited  amount  of  research  into  naval 
applications  of  MOS.  In  particular,  statistical  studies  have 
been  done  into  forecasting  Levante  winds  in  Spain  [Godfrey  and 


Lowe,  1979],  ceiling  and  visibility  prediction  in  the  southern 
California  (SOCAL)  naval  operating  area  [Lewit,  1980],  marine 
fog  and  visibility  predictability  in  the  North  Pacific  Ocean 
[Renard  et  al.,  1983  and  Renard  and  Thompson,  1984],  Presently, 
a  program  is  in  operation  which  provides  MOS  forecasts  for 
selected  U.S.  Navy  and  Marine  Corps  CONUS  locations.  These 
services,  which  are  made  available  from  NWS,  are  based  on  the 
National  Meteorological  Center  (NMC)  limited  fine  mesh  model 
predictions.  This  MOS  program  was  initiated  on  10  November 
1982  and  provides  forecasts  for  twelve  weather  parameters 
which  include  visibility,  obstructions  to  visibility  and 
cloud  amount  [Naval  Environmental  Prediction  Research 
Facility,  1982] . 

The  results  of  these  limited  studies  along  with  the 
encouraging  performances  of  both  the  NWS  and  AWS  MOS  programs 
and  the  implementation  of  the  Navy  Operational  Global  Atmos¬ 
pheric  Prediction  System  (NOGAPS)  dynamical  primitive  equation 
(PE)  model  at  the  Fleet  Numerical  Oceanography  Center  (FNOC) , 
in  Monterey,  California  prompted  the  decision  in  September 
1982  for  the  Navy  to  pursue  its  own  MOS  program. 

Fig.  1  is  an  overview  of  the  currently  proposed  milestones 
for  the  Navy  MOS  program.  The  first  operational  weather 
parameter  being  investigated  in  this  proposed  ten-year  Navy 
effort  is  horizontal  visibility  at  sea,  with  the  initial  goal 
of  this  project  being  the  investigation  and  development  of 
statistical  predictive  schemes  for  forecasting  horizontal 
visibility  over  the  North  Atlantic  Ocean. 


The  impact  of  fog  and  other  impediments  to  visibility  on 
naval  operations  is  well  documented  throughout  maritime 
history.  Records  show  countless  catastrophes  and  accidents 
which  were  directly  attributable  to  poor  visibility  at  sea. 

For  example,  on  29  May  1914,  the  Canadian  liner  Empress  of 
Ireland  collided  with  the  Norwegian  vessel  Storstad  in  dense 
fog  on  the  Saint  Lawrence  River  resulting  in  1,024  fatalities 
and  similarly,  the  legendary  "North  Sea  haze"  was  a  critical 
element  in  the  World  War  I  tactics  employed  at  Jutland  in 
1916.  Also,  one  of  the  most  spectacular  maritime  disasters 
in  the  U.S.  Navy's  history  took  place  on  9  September  1923 
when  seven  Pacific  fleet  destroyers  struck  the  rocks  and  ran 
aground  in  dense  fog  off  of  Point  Arguello,  California. 

Research  into  predicting  marine  visibility  via  traditional 
linear  regression  methodologies  has  taken  place  at  the  Naval 
Postgraduate  School  (NPS)  since  the  early  1960's.  Generally, 
early  visibility  forecasting  experiments  identified  potential 
physical  air/ocean  mechanisms  [Schramn,  1966]  and  emphasized 
the  inherent  likelihood  of  human  error  in  at-sea  visibility 
observations  [Nelson,  1972] .  Later  experimentation  by 
Aldinger  (1979),  Yavorsky  (1980)  and  Selsor  (1980)  concen¬ 
trated  on  various  modifications  to  multiple  linear  regression 
schemes  and  the  analysis  of  prediction  skill  measurements. 

This  study  presents  a  direct  follow-on  to  the  research 
presented  by  Karl  (1984),  in  which  statistical  methodologies 
presented  by  Preisendorf er  (1983a, b,c)  and  multiple  linear 


regression  techniques  presented  by  Lowe  (1984a)  were  compared 
and  contrasted.  In  Karl's  preliminary  study,  Preisendorf er ' s 
three  strategies,  two  based  on  maximum  conditional  proba¬ 
bility  and  one  based  on  natural  regression,  as  well  as 
Lowe's  linear  regression  threshold  models  were  tested  and 
applied  to  sets  of  FNOC  model  output  parameters  (MOPs)  from 
both  the  North  Pacific  and  North  Atlantic  Ocean  areas.  The 
North  Atlantic  Ocean  study  was  separated  into  effective 
physically  homogeneous  areas  [Lowe,  1984b].  Karl's  study 
specifically  dealt  with  an  evaluation  of  the  MOS  scheme 
applied  to  oceanic  regions  for  the  TAU-00  model  output 
during  the  period  15  May  to  07  July  1983. 

This  study  concerns  itself  with  a  continued  evaluation 
and  further  refinement  of  statistical  methods  proposed  by 
Preisendorfer  as  well  as  the  linear  regression  threshold 
models  presented  by  Lowe.  With  reference  to  Karl's  study, 
other  North  Atlantic  Ocean  areas  and  model  forecast  projections 
(e.g.  TAU-24  and  TAU-48)  are  addressed. 


OBJECTIVES  AND  APPROACH 


The  primary  objectives  of  this  study  are  to  continue  the 
previous  NPS  horizontal  marine  visibility  prediction  re¬ 
search  initiated  by  Karl  (1984)  and  to  continue  the  search 
for  an  optimal  Model  Output  Statistics  (MOS)  prediction 
scheme  to  operationally  forecast  coastal  and  open  ocean 
visibility  over  the  North  Atlantic  Ocean.  The  approach 
employed  in  meeting  the  stated  objectives  is  listed  below: 

A.  Apply  and  evaluate  the  Preisendorf er  maximum  proba¬ 
bility  and  natural  regression  strategies  (1983a, b,c)  to  addi 
tional  North  Atlantic  Ocean  homogeneous  areas  [Lowe,  1983b] 
using  May  through  July  1983,  NOGAPS  predictand/predictor 
data . 

B.  Expand  the  Model  Output  Predictor  (MOP)  data  sets  to 
include  the  NOGAPS  model  TAU-00,  and  the  TAU-24  and  TAU-48 
prognostic  times  defined  in  Chapter  III. 

C.  Investigate  specific  two-stage,  equal  variance  and 
quadratic  multiple  linear  regression  threshold  models  pro¬ 
posed  by  Lowe  (1984a)  for  the  oceanic  areas  and  model  output 
periods  addressed  in  A.  and  B.  above. 

D.  Compare  and  contrast  the  individual  results  of  the 
Preisendorf er  statistical  methodologies  to  those  of  the  Lowe 
approach . 


E.  Conduct  a  limited  series  of  experiments  in  which  a 


1984  data  set,  15  May  to  23  June,  is  utilized  as  an 


independence  evaluation  of  the  predictive  models  constructed 
with  1983  NOGAPS  data. 

F.  Based  on  A.  to  E.  above,  present  an  interim  recommenda¬ 
tion  for  an  optimal  statistical  approach  to  forecast  North 
Atlantic  Ocean  horizontal  visibility  as  a  function  of  pre¬ 
diction  time  and  homogeneous  area. 


III.  DATA 


A.  VISIBILITY  OBSERVATIONS  AND  SYNOPTIC  CODES 

Horizontal  visibility  observations  taken  from  seagoing 
platforms  are  reported  as  values  of  ten  standardized  World 
Meteorological  Organization  (WMO)  synoptic  weather  codes. 
These  codes  range  in  value  from  90,  which  corresponds  to 
visibility  less  than  50  meters,  to  99,  which  corresponds  to 
visibility  equal  to  or  greater  than  50  kilometers.  Human 
observational  error  and  inexactness  in  measuring  visibility 
at  sea  necessitates  a  generalization  of  visibility  classifi¬ 
cation  for  prediction  purposes,  as  follows: 

Visibility  Category  Synoptic  Code  Visibility  Range 


I 

90-94 

<  2  km 

II 

95-96 

>  2  km  to  <  10  km 

III 

97-99 

>  10  km 

The  above  scheme  coincides  with  the  classification  scheme 
proposed  by  Karl  (1984)  and  is  based  upon  the  below  listed 
U.S.  Navy  operational  criteria. 

1.  10  km  (5  n  mi) — U.S.  Navy  aircraft  carrier  at-sea 
flight  recovery  operations  change  from  visual  (VFR)  to 
controlled  (IFR)  approach  guidelines  [Department  of  the 
Navy ,  1979] . 

2.  2  km  (1  n  mi) --the  sounding  of  reduced  visibility 
signals  for  all  vessels  operating  in  international  waters. 


The  term  "reduced  visibility"  is  not  specifically  defined  in 
the  International  Regulations  for  Preventing  Collisions  at 
Sea,  1972.  The  distance  of  1  n  mi  is  generally  considered 
to  be  the  governing  operational  distance. 

B.  NORTH  ATLANTIC  OCEAN  DATA 

1 .  Area 

The  North  Atlantic  Ocean,  from  0°-80°  N  latitude, 
was  divided  into  homogeneous  oceanic  areas  by  Lowe  (1984b) 
using  a  statistical  cluster  analysis  technique.  The  specific 
homogeneous  areas  evaluated  in  this  study  are  identified  as 
areas  2,  3W  and  4  on  Fig.  2.  These  areas  were  selected  be¬ 
cause  they  individually  represent  a  range  of  different  rela¬ 
tive  frequencies  of  poor  visibility  observations.  Area  3W, 
which  was  used  by  Karl  (1984)  for  his  preliminary  experimen¬ 
tation,  represents  an  area  of  relatively  frequent  occurrence 
of  poor  visibility,  while  area  4  represents  an  area  of  rela¬ 
tively  sparse  occurrence  of  poor  visibility  and  area  2 
represents  an  intermediate  case. 

2 .  Time  Period 

Data  from  mid-May  1983  to  mid-July  1983  were  combined 
to  form  a  more  extensive  data  set,  hereafter  referred  to  as 
FATJUNE  1983.  FATJUNE  1983  was  selected  as  the  initial  data 
set  for  statistical  experimentation  because  of  the  high  fre¬ 
quency  of  occurrence  of  poor  visibility  observations  for 
many  areas  of  the  North  Atlantic  Ocean  during  this  period. 


1200  GMT  synoptic  ship  report  data  were  used  exclusively  in 
this  study.  This  time  corresponds  to  general  daylight  condi¬ 
tions  over  the  North  Atlantic  Ocean  during  FATJUNE.  In 
addition  to  FATJUNE  1983,  a  limited  May  15  to  June  23  1984 
data  set,  possessing  the  same  geographical  coverage  and  day¬ 
light  characteristics  of  FATJUNE  1983,  was  utilized  in  an 
independent  test  of  the  predictive  arrays  and  equations 
generated  in  this  study. 

For  the  purpose  of  this  study,  TAU-00  generally 
represents  six-hour  model  forecast  fields  valid  at  1200  GMT. 
Three  specific  fields,  namely  temperature,  geopotential 
height  and  wind,  are  model  initialization  fields  valid  at 
1200  GMT.  TAU-24  and  TAU-48  are  defined  as  24-hour  and 
48-hour  model  forecast  fields,  valid  at  1200  GMT.  TAU-00, 
TAU-24  and  TAU-48  model  output  parameters  (predictors)  are 
employed  in  the  00,  24  and  48  hour  forecast  schemes,  respect¬ 
fully.  Summaries  of  the  visibility  frequencies  for  each 
visibility  category,  as  a  function  of  homogeneous  area 
and  prediction  time,  for  FATJUNE  1983  and  the  15  May  to  23 
June  1984  data  set,  are  contained  in  Tables  I  through  III 
and  Table  VI  respectively. 

3 .  Synoptic  Weather  Reports 

All  synoptic  visibility  observations  (predictand  data) 
for  this  study  were  provided  by  the  Naval  Oceanography  Com¬ 
mand  Detachment  (NOCD) ,  Asheville,  North  Carolina  which  is 
co-located  with  the  National  Climatic  Data  Center  (NCDC) . 
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The  observations  which  contained  systematic  observer  error 
or  were  obviously  erroneous,  as  determined  from  the  data 
quality  indicators  provided  with  the  data,  were  deleted 
from  the  working  data  sets. 

4 .  Predictor  Parameters 

Fifty  TAU-00,  fifty-four  TAU-24  and  fifty-four  TAU-48 
model  output  predictors  (MOP's)  were  provided  by  the  Fleet 
Numerical  Oceanography  Center  (FNOC) ,  Monterey,  California. 
These  parameters  are  generated  by  their  current  operational 
atmospheric  prediction  model,  the  Navy  Operational  Global 
Atmospheric  Prediction  System  (NOGAPS).  All  MOP’s  were 
interpolated  from  model  grid  coordinates  to  synoptic  ship 
report  position  using  a  linear  interpolation  scheme.  In 
addition  to  the  initial  group  of  model  output  parameters, 
ten  derived  parameters  representing  calculated  quantities, 
such  as  parameter  gradients  and  products,  were  included  as 
potential  predictors.  Of  the  entire  group  of  potential 
predictor  parameters,  only  forty  TAU-00  and  forty-seven 
TAU-24  and  TAU-48  MOP's  were  actually  used  to  develop  the 
various  Preisendorfer  (  1983a, b,c)  and  linear  regression 
threshold  models  [Lowe,  1983a].  The  remainder  of  the  NOGAPS 
model  output  parameters  were  dropped  from  consideration  because 
1)  the  MOP  lacked  a  physical  linkage  to  the  visibility  pre- 
dictand  and/or  2)  a  lack  of  significant  digits  (lost  during 
the  transfer  of  FNOC  data  to  the  main  computer  center's 
mass  storage  system)  rendered  the  particular  MOP  useless. 
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A  list  of  all  available  TAU-00 ,  TAU-24  and  TAU-48  MOP's  are 


> 


included  in  Appendix  B. 

For  each  homogeneous  area  and  model  forecast  projec¬ 
tion,  a  set  of  three  linear  regression  equations,  in  addition 
to  the  aforementioned  MOP's,  were  included  as  potential 
MOP's  for  a  separate  evaluation  of  the  Preisendorf er  methodology 
(the  PR+-BMD  model).  These  three  predictor  equations  were 
obtained  from  two  standardized  linear  regression  software 
packages,  namely  P2R — stepwise  regression  and  P9R — all 
possible  subsets  regression,  as  addressed  in  the  BMDP  Sta¬ 
tistical  Software  [University  of  California,  1983].  The 
P2R  was  initially  employed  in  the  evaluation  of  areas  2  and 
4,  TAU-00  data,  while  the  P9R  program  was  employed  in  the 
remainder  of  the  cases  studied.  The  change  to  the  P9R 
program  was  initiated  as  a  safeguard  against  any  potential 
predictor  selection  bias  incorporated  in  the  P2R  software. 
Specific  details  concerning  these  statistical  software 
packages  are  addressed  in  Appendix  A. 

C.  DEPENDENT/INDEPENDENT  DATA  SETS 

Due  to  the  limited  amount  of  data  available  to  this 
study  for  each  of  the  North  Atlantic  Ocean  homogeneous  areas, 
it  was  necessary  to  withhold  a  significant  amount  of  the 
observations  from  the  developmental  model  to  use  as  an 
independent  data  set.  That  amount  was  set  as  one-third  for 
the  experiments  reported  here.  This  was  accomplished  by  the 
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use  of  a  counter  and  transfer  statement  in  the  computer 
programs  which  prevented  every  third  observation  from  enter¬ 
ing  the  developmental  computations.  To  ensure  that  the 
dependent  and  independent  data  were  representative  of  the 
same  population,  a  95%  confidence  interval  for  proportions 
[Miller  and  Freund,  1977]  was  established  from  the  entire 
data  set,  for  each  visibility  category,  and  the  dependent 
and  independent  data  sets  were  constrained  to  have  visibility 
frequencies  within  these  established  confidence  intervals. 
Table  IV  summarizes  the  dependent  and  independent  data  for 
the  North  Atlantic  Ocean  data  set. 
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IV.  PROCEDURES 


A.  TERMS  AND  SYMBOLS 

The  terms  and  statistical  symbols  defined  below  will 
be  used  throughout  the  remainder  of  this  report.  The 
formal  mathematical  definitions  are  described  in  Karl  (1984) 

1.  Maximum  probability  strategy — choosing  forecast 
visibility  category  based  upon  the  highest  conditional 
probability  of  visibility  within  a  predictor  interval. 

a.  MAXPROB  I — designation  of  the  maximum  probability 
strategy  in  which  ties  of  the  highest  conditional  probabili¬ 
ties  in  a  predictor  interval  are  resolved  by  the  generation 
of  a  random  number. 

b.  MAXPROB  II — designation  of  the  maximum  probability 
strategy  in  which  ties  of  the  highest  conditional  probabili¬ 
ties  in  a  predictor  interval  are  resolved  by  assigning  the 
lowest  visibility  category,  of  those  tied,  as  the  forecast 
category . 

2.  Natural  regression  strategy-choosing  forecast  visi¬ 
bility  categories  based  upon  the  statistical  average  of  the 
conditional  probabilities  of  visibility  within  a  predictor 
interval . 

3.  AO — the  probability  of  a  zero-class  visibility  cate¬ 
gory  forecast  error  (e.g.,  if  visiblity  category  I  is 
forecast  and  observed) . 


4.  Al — the  probability  of  a  one-class  visibility  category 
forecast  error  (e.g.,  if  visibility  category  I  is  forecast 
and  category  II  is  observed)  . 

5.  CE — class  error  parameter  defined  as  A0+2A1,  used  as 
the  primary  aid  in  identifying  the  first  predictor. 

6.  PP — the  potential  predictability  of  visibility  by  any 
given  predictor. 

7.  Functional  dependence.  This  is  a  measure  of  the 
stochastic  dependence  of  one  predictor  upon  another.  Func¬ 
tional  dependence  is  the  probability  that  one  of  the  predic¬ 
tors  will  change  when  the  other  does.  High  functional 
dependence  values  between  one  already  selected  predictor  and 
another  potential  predictor,  indicates  that  little  addi¬ 
tional  information  beyond  the  selected  predictor  is  possible. 
Conversely,  a  low  functional  dependence  value  between  the 
same  two  p vad  "'.ctors ,  indicates  that  each  predictor  possesses 
a  high  degree  of  linearly  uncorrelated  information  concerning 
the  predictand.  Functional  dependence  range  is  0.0  to  1.0 
(1.0  =  highest  functional  dependence).  The  specific  deriva¬ 
tion  and  mathematical  description  of  the  concept  of  "func¬ 
tional  dependence"  is  discussed  in  greater  depth  by 
Preisendorfer  (1983c). 

8.  Root-sum- squared  functional  dependence.  The  functional 
dependence  of  a  predictor  on  all  predictors  already  included 
in  the  developmental  model.  It  is  equal  to  the  square-root 


of  the  sum  of  the  squares  of  the  individual  functional 
dependence  values. 

9.  TS1 — threat  score  for  visibility  category  I,  as 
computed  from  a  contingency  table  (see  Appendix  C) . 

10.  ATS1 — adjusted  threat  score  for  visibility  category 
I  which  removes  the  influence  of  the  data  set  category 
frequency  (see  Appendix  C) . 

11.  AAO — adjusted  AO.  A  contingency  table  statistic 
which  removes  the  influence  of  the  most  frequent  visibility 
category  in  a  set  of  data  (similar  to  a  normalized  value) 
(see  Appendix  C) . 

B .  COMPUTER  PROGRAMS 

Four  computer  programs  were  developed  to  test  the  pro¬ 
posed  Preisendorfer  (1983a, b,c)  methodology.  The  programs 
are  on  file  in  the  Department  of  Meteorology,  Naval  Post¬ 
graduate  School,  Monterey,  California,  93943. 

1.  A  program  to  compute  AO,  Al,  CE  and  PP  for  all  predic¬ 
tors,  all  strategies  (MAXPROB  I,  MAXPROB  II  and  natural 
regression)  for  a  particular  number  of  equally  populous 
predictor  intervals.  Statistics  for  the  three  strategies 
are  based  upon  the  same  predictor (s)  rather  than  the  best 
predictor (s)  for  each  strategy. 

2.  A  program  to  compute  functional  dependence  values  for 
all  predictors,  on  a  given  predictor,  for  a  given  number  of 
equally  populous  predictor  intervals  and  to  compute  the 


associated  96%  critical  confidence  interval  value,  referred 
to  as  functional  dependence (9 6)  in  this  study,  by  Monte  Carlo 
means . 

3.  A  program  to  construct  contingency  tables  and  to 
compute  skill  and  threat  scores,  for  both  the  dependent  and 
independent  data  sets . 

4.  A  program  to  generate  100  random  data  sets,  from  the 
marginal  probabilities  of  the  predictor (s)  in  the  develop¬ 
mental  model,  and  to  compute  upper  and  lower  5%  critical 
confidence  interval  values  for  AO  and  Al  to  be  used  for 
testing  the  significance  of  the  results  from  each  of  the 
Preisendorfer  models  against  chance.  These  confidence 
interval  values  are  calculated  via  Monte  Carlo  means. 

C .  MODELS 

1 .  Preisendorfer  PR  Model 

This  model  represents  the  first  of  two  different 
applications  of  the  basic  Preisendorfer  methodology 
[Preisendorfer,  1983a, b,c].  Karl  (1984),  in  his  preliminary 
research,  provides  a  rigorous  interpretation  and  results 
associated  with  this  statistical  forecasting  methodology. 
Karl's  study  provides  the  necessary  background  for  the  con¬ 
tinued  investigation  and  evaluation  of  this  model  and  readers 
interested  in  specific  details  are  advised  to  consult  this 
document . 

The  PR  model  utilizes  the  working  set  of  NOGAPS 
model  output  parameters  (MOP's)  and  derived  parameters 


(Appendix  B)  as  potential  predictors  in  constructing  a 
developmental  model,  based  upon  the  dependent  data  set, 
which  provides  the  structure  by  which  the  independent  data 
set  is  tested  and  evaluated.  In  general,  these  potential 
predictors  have  their  range  of  values  partitioned  into 
discretized  equally  populous  predictor  intervals  ("cells") 
and  conditional  probabilities  of  the  predictand  are  calcu¬ 
lated  according  to  the  three  modified  visibility  categories 
(VISCAT)  I,  II  and  III.  Three  separate  strategies  of  deter¬ 
mining  the  specific  VISCAT  which  is  to  be  identified  with 
each  predictor  value,  are  proposed.  These  strategies,  two 
based  upon  maximum  probability  and  the  third  based  on  a 
natural  regression  approach,  are  addressed  as  MAXPROB  I, 
MAXPROB  II  and  natural  regression  in  the  remaining  portions 
of  this  study. 

Initial  evaluation  of  this  model  involves  varying  the 
equally  populous  predictor  intervals  from  sizes  of  four 
through  ten,  and  selecting  an  optimal  first  predictor  which 
provides  one  of  the  following  requirements  in  the  designated 
order : 

a.  the  lowest  CE  value  of  all  the  potential  predictors 

b.  the  highest  PP  value  of  all  the  potential  predictors 

Once  a  first  predictor  is  identified  for  each  of  the 
four  through  ten  equally  populous  predictor  intervals, 
corresponding  VISCAT  I,  II  and  III  threat  and  AO  skill 
scores  (Appendix  E)  are  calculated  for  both  the  dependent 


and  independent  data  sets.  The  practice  of  selecting  an 
optimal  equally  populous  predictor  interval  from  the  eligible 
grouping  sizes  of  four  through  ten,  was  proposed  by  Karl 
(1984)  as  a  practical  procedure  which  would  permit  the 
realization  of  peak  skill  scores  as  well  as  maintain  asso¬ 
ciated  computer  storage  requirements  at  a  manageable  level. 

An  unfortunate  consequence  of  this  range  of  potential  group¬ 
ing  sizes  is  that  certain  statistical  calculations  associated 
with  equally  populous  predictor  intervals  of  eight,  nine 
and  ten  are  terminated  before  completion  due  to  a  two  mega¬ 
byte  storage  ceiling  at  the  NPS  W.R.  Church  Computer  Center. 
When  considering  potential  predictor  intervals,  the  size  of 
the  interval  is  of  obvious  importance,  with  lower  values 
being  the  most  desirable.  The  criterion  for  determining  the 
optimal  equally  populous  predictor  interval  is  to  select  the 
smallest  interval  value  which  maximizes  the  dependent  data 
set  adjusted  AO  and  independent  adjusted  VISCAT  I  threat 
score.  For  this  study,  this  interval  value  was  fixed  for 
all  ensuing  aspects  of  the  model  evaluation.  In  practice, 
the  selection  of  equally  populous  predictor  intervals  was 
based  upon  the  initial  adjusted  AO  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
MAXPROB  II  strategy.  The  MAXPROB  II  scores  were  routinely 
found  to  be  the  highest  for  each  case  evaluated,  at  this 
early  stage  in  the  evaluation  process,  and  therefore  used 
as  the  basis  for  grouping  selection.  As  the  equally  populous 


grouping  interval  remains  constant  throughout  the  Preisen- 
dorfer  models,  the  MAXPROB  I  and  natural  regression  strate¬ 
gies  practically  play  no  role  in  the  predictor  selection 
process . 

Once  the  first  predictor  and  its  associated  equally 
populous  predictor  interval  have  been  identified,  a  func¬ 
tional  dependence  test  of  the  first  predictor  against  those 
remaining  potential  predictors  is  run.  The  second,  third 
and  all  subsequent  predictors  are  selected  only  if  both  of 
the  following  criteria  are  met: 

a.  subsequent  predictors  must  increase  AO  over  the 
AO  value  attained  at  the  preceding  level,  and 

b.  the  selected  predictor  must  have  the  lowest 
functional  dependence  and  root-sum- square 
functional  dependence  of  all  the  remaining 
potential  predictors. 

After  each  predictor  selection  stage  has  been  com¬ 
pleted,  significance  tests  are  run  upon  the  developmental 
model  to  determine  if  the  results  are  suitably  significant 
as  compared  to  random  chance.  This  testing  is  accomplished 
via  Monte  Carlo  testing  methods  using  the  conditional 
probabilities  of  the  selected  predictors  and  assuming  equal 
probability  of  occurrence  for  the  three  modified  visibility 
categories.  Functional  dependence/root-sum-square  functional 
dependence,  AO,  and  A1  statistics  are  calculated  for  each  of 
100  randomly  generated  data  sets.  For  the  developmental 
model  to  yield  results  which  are  significant  at  the  speci¬ 
fied  confidence  interval  values,  each  one  of  the  following 
criteria  must  be  met: 


a.  AO  must  be  equal  to  or  greater  than  AO (96) 

b.  Al  must  be  equal  to  or  less  than  Al(05) 

c.  the  functional  dependence  value  for  a  selected 
predictor  must  be  less  than  functional 
dependence (96) 

As  with  the  process  of  selecting  equally  populous 
predictor  intervals,  the  AO,  A0(96),  Al  and  Al(05)  statistics 
(Appendix  G) ,  reflect  scores  for  the  MAXPROB  II  strategy. 

The  AO  statistics  routinely  were  found  to  be  the  highest  for 
thif.  trategy  and  thus  were  used  as  the  basis  for  ensuring 
the  aforementioned  predictor  selection  criteria  were  met. 
However,  the  MAXPROB  I  strategy  often  produced  AO  values 
identical  to  MAXPROB  II.  The  natural  regression  strategy 
regularly  lagged  the  two  maximum  probability  strategies  in 
AO  and  Al  scores  and  consequently  played  no  real  role  in  the 
prediction  selection  process.  Specific  trends  in  A0/A1 
scores  can  be  seen  in  Appendix  G. 

From  a  practical  standpoint,  the  model  development 
continues  until  computer  storage  limitations  preclude  further 
addition  of  predictors.  This  generally  occurred  at  the  fifth 
predictor  level. 

Once  the  developmental  model  is  completed,  contingency 
tables  of  forecast  visibility  category  versus  observed  visi¬ 
bility  category  are  constructed  for  both  the  dependent  and 
independent  data  sets,  and  threat  and  skill  scores  are 
computed  and  compared. 


2. 


Preisendorfer  PR+BMD  Model 


This  model  is  still  the  PR  model  described  above. 

Now,  sets  of  three  linear  regression  equations  (Appendix  D) 
are  added  to  the  list  of  potential  NOGAPS  model  output  and 
derived  predictor  parameters.  The  inherent  difference  of 
these  predictors  is  evidenced  in  both  the  predictor  selec¬ 
tion  process  as  well  as  in  the  resulting  skill  and  threat 
scores,  as  will  be  demonstrated  in  Chapter  V. 

3 .  Equal  Variance  Threshold  Model  (EVAR) 

This  model  represents  the  first  of  two  threshold 
models,  developed  by  Lowe  (1984a),  which  were  evaluated  in 
this  study.  The  model  uses  an  algorithm  which  requires  the 
assumption  that  the  variances  of  two  normally  distributed 
populations  which  are  to  be  separated  by  a  threshold  value 
are  equal,  while  their  means  are  unequal.  A  detailed  dis¬ 
cussion  of  the  theoretical  background  of  this  scheme  is 
addressed  in  Appendix  A. 

A  two-stage  separation  scheme  was  used  to  effectively 
divide  the  visibility  categories  (VISCAT)  I,  II  and  III 
into  a  first-stage  VISCAT  I  versus  a  combined  VISCAT  II  plus 
VISCAT  III  separation,  and  subsequently  VISCAT  II  versus 
VISCAT  III  separation  for  each  homogeneous  area  and  model 
output  time.  This  separation  was  accomplished  by  setting 
all  VISCAT  I  observations  equal  to  an  arbitrary  integer  value 
of  zero  and  the  combined  VISCAT  II  plus  VISCAT  III  observations 
equal  to  an  arbitrary  integer  value  of  one  and  generating 


a  linear  regression  equation  to  suitably  describe  the 
resulting  two  distributions.  This  linear  regression  equation 
was  then  used  in  the  graphical  plotting  program  BMDP5D, 
from  the  BMDP  Statistical  Software  [University  of  California, 
1983] ,  to  generate  a  set  of  histograms  describing  the  first 
stage  separation.  Included  with  the  graphical  histogram 
output  is  a  listing  of  the  individual  frequency  of  observa¬ 
tion  (P) ,  mean  (y)  and  standard  deviation  (a)  of  each  of  the 
specified  visibility  distributions.  These  statistics  are 
incorporated  into  the  equal  variance  threshold  algorithm 
and  a  corresponding  threshold  value  is  calculated. 

Following  the  first-stage  threshold  calculation,  a 
second  linear  regression  equation  is  generated,  based  upon 
only  those  VISCAT  II  plus  VISCAT  III  observations  which 
exceed  the  previously  calculated  threshold  value.  This 
effectively  eliminates  any  VISCAT  II  plus  VISCAT  III  obser¬ 
vations  less  than  the  threshold  value  (i.e.,  those  observa¬ 
tions  contained  in  the  tail  of  the  distribution) ,  from  being 
included  in  the  second-stage  regression.  The  previous  proce¬ 
dure  of  generating  corresponding  histograms  and  statistics 
is  repeated,  based  upon  all  VISCAT  II  observations  being 
assigned  an  arbitrary  integer  value  of  zero  and  all  VISCAT 
III  observations  being  assigned  an  integer  value  of  one.  A 
second-stage  equal  variance  threshold  value  is  then  calcu¬ 
lated  which  separates  VISCAT  II  from  VISCAT  III. 
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With  the  two-stage  separation  complete,  the  indepen¬ 
dent  data  set  is  processed  through  the  governing  equations 
and  thresholds  to  obtain  a  set  of  observed  visibility  value 
results  versus  calculated  "forecast"  visibility  value  re¬ 
sults.  These  results,  in  contingency  table  format  for  each 
evaluated  case,  are  presented  in  Chapter  V  and  Appendix  G. 

4 .  Quadratic  Threshold  Model  (QUAD) 

This  model  represents  the  second  of  two  threshold 
models,  developed  by  Lowe  (1984a),  which  were  evaluated  in 
this  study.  The  model  uses  an  algorithm  which  requires  the 
assumption  that  both  the  variances  and  the  means  of  two 
normally  distributed  populations,  which  are  to  be  separated 
by  a  threshold  value,  are  equal.  Similar  to  the  EVAR  model, 
a  detailed  discussion  of  the  theoretical  background  of  this 
scheme  is  addressed  in  Appendix  A. 

The  general  two-stage  separation  procedure  employed 
with  this  model  is  identical  to  that  described  for  the  EVAR 
model  in  IV. C.  above.  The  only  difference  between  the  QUAD 
and  EVAR  model  is  the  algorithm,  based  upon  a  solution  to 
a  quadratic  equation  in  this  model,  used  to  calculate  the 
appropriate  threshold  values. 

5 .  Maximum-Likelihood-of- Detect ion  Model 

The  maximum-likelihood-of-detection  criteria  (MLDC) 
is  an  additional  threshold  technique  which  is  included  in 
this  study  as  a  possible  alternative  to  the  aforementioned 
EVAR  and  QUAD  minimum  probable  error  threshold  models.  The 
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MLDC  involves  calculating  a  threshold  value  based  upon  the 
assumptions  that  the  population  frequencies  and  variances 
of  two  normally  distributed  samples  which  are  to  be  separated 
are  identical.  This  technique  is  particularly  well  suited 
for  cases  where  the  threat  frequency  (i.e.,  number  of 
threatening  events  divided  by  the  total  number  of  threat  and 
non-threat  events)  approaches  very  small  values  (e.g., 
statistical  rare  events) . 

Unlike  the  EVAR  and  QUAD  models,  the  two-stage 
separation  employed  with  this  technique  utilizes  a  first- 
stage  VISCAT  I+II  versus  VISCAT  III  followed  by  a  second- 
stage  VISCAT  I  versus  VISCAT  II  separation.  In  calculating 
the  specific  threshold  values,  the  lowest  frequency  visibility 
category  (usually  the  VISCAT  I  threat  category)  is  assigned 
an  arbitrary  integer  value  of  one.  The  remaining  larger 
visibility  category/ies  are  assigned  the  arbitrary  integer 
value  of  zero.  Proceeding  in  the  same  manner  as  described 
with  the  EVAR  and  QUAD  models,  population  means  are  calcu¬ 
lated  for  each  separation  stage.  The  threshold  value  is 
simply  the  mid-point  between  the  two  population  means.  A 
detailed  discussion  of  the  theoretical  background  of  this 
scheme  is  addressed  in  Appendix  A. 


V.  RESULTS 


The  general  procedures  outlined  in  Chapter  I  were  fol¬ 
lowed  in  evaluating  the  statistical  scoring  techniques  for 
the  oceanic  homogeneous  areas  2,  3W  and  4.  Certain  slight 
modifications  were  required  to  handle  the  relatively  low 
frequency  of  visibility  category  I,  in  area  4  for  the  TAU- 
00,  TAU-24  and  TAU-48  model  output  data  sets.  Fig.  2 
displays  the  individual  oceanic  homogeneous  areas  for  FATJUNE 
1983.  Tables  I  through  III  identify  the  frequency  of  occur¬ 
rence  of  visibility  categories  I,  II  and  III  at  TAU-00, 

TAU-24  and  TAU-48  for  each  of  the  evaluated  homogeneous 
areas . 

In  discussing  the  results  of  this  study,  specific  comment 
is  focused  upon  the  optimal  model  for  each  case  as  well  as 
any  significant  finding  observed  by  the  author.  Certain 
characteristics  of  the  evaluated  cases  are  repetitious  and 
are  considered  adequately  described  by  their  associated 
figures.  Consequently,  the  entire  assemblage  of  figures  in 
Appendix  G  are  not  individually  addressed.  These  figures 
are  nevertheless  considered  noteworthy,  as  they  document  the 
performance  of  each  tested  model  in  this  study,  and  are 
included  as  a  matter  of  record.  The  following  presentation 
of  the  results  of  this  experimentation  are  arranged  accord¬ 
ing  to  the  specific  oceanic  homogeneous  area  and  model 
output  period. 


In  general,  four  models  are  evaluated  for  each  of  the 
predefined  homogeneous  areas/model  forecast  projections.  The 
four  models  are:  the  Preisendorf er  methodology  utilizing 
NOGAPS  model  output  predictors  and  a  limited  number  of 
derived  predictors  (PR) ,  the  Preisendorf er  methodology  uti¬ 
lizing  both  NOGAPS  model  output  predictors,  derived  predic¬ 
tors  and  linear  regression  equation  predictors  (PR+BMD),  an 
equal  variance  linear  regression  threshold  model  (EVAR)  and 
a  quadratic  linear  regression  threshold  model  (QUAD) . 

A.  NORTH  ATLANTIC  OCEAN,  AREA  2 

Area  2  encompasses  a  geographic  region  that  extends  from 
the  southeastern  tip  of  Newfoundland,  across  the  North 
Atlantic  Ocean  to  the  eastern  coast  of  England,  north  to 
the  Five  Fingers  of  Iceland  and  back  to  the  Canadian  coast 
north  of  Newfoundland.  Fig.  2  gives  the  pictorial  repre¬ 
sentation  of  the  area. 

1 •  Area  2,  TAU-00 

Fig.  3  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  £0  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR  model.  For  this  case,  a  grouping  size  of  eight  was  se¬ 
lected.  Results  of  the  individual  MAXPROB  I,  MAXPROB  II 
and  Natural  Regression  strategies  are  shown  in  Figs.  4a 
though  4c.  The  MAXPROB  II  strategy  (Fig.  4b)  produced  the 
largest  overall  independent  data  VISCAT  I  adjusted  threat 
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score,  namely  0.23  (unadjusted,  0.30).  This  peak  threat 
score  occurs  with  the  inclusion  of  the  first  predictor, 

E850,  and  declines  marginally  with  the  addition  of  the  re¬ 
maining  four  predictors.  Of  the  three  strategies,  the 
natural  regression  strategy  (Fig.  4c),  yields  the  poorest 
overall  threat  scores  with  its  peak  threat  score  occurring 
with  the  addition  of  the  fourth  predictor.  The  predictors 
selected  for  this  case  are  E850,  ENTR,  DVDP,  U1000,  and 
STRTH . 

The  associated  functional  dependence  and  A0/A1  sta¬ 
tistics  and  96%/05%  confidence  interval  values  for  these 
predictors  are  shown  in  Fig.  5.  The  trend  of  functional 
dependence  versus  its  96%  confidence  interval  shows  that  the 
specific  functional  dependence  values  associated  with  the 
chosen  predictors  never  falls  within  the  96%  confidence 
interval.  At  the  first  predictor  level,  for  example,  the 
functional  dependence  of  ENTR  upon  E850  has  a  value  of  0.1146 
as  compared  to  a  96%  confidence  interval  value  of  0.1039. 

This  infers  that  the  corresponding  scores  (i.e.,  threat 
scores,  A0  and  Al)  are  not  statistically  significant  at  the 
preselected  96%  confidence  interval  level. 

Fig.  6  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  A0  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR+BMD  model.  For  this  case  an  equally  populous  grouping 
size  of  seven  was  selected.  Results  of  the  three  individual 
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Preisendorf er  strategies,  along  with  the  corresponding  con¬ 
tingency  tables,  can  be  seen  in  Figs.  7a  through  7c.  As 
with  the  PR  model,  a  maximum  independent  VISCAT  I  threat 
score  was  obtained  with  the  MAXPROB  II  strategy  using  the 
first  predictor  selected,  namely  the  linear  regression  equa¬ 
tion  predictor  BMDl  (Appendix  D) .  The  overall  independent 
adjusted  VISCAT  I  threat  score  achieved  with  this  model  is 
0.29  (unadjusted,  0.36),  which  is  .06  greater  than  that  for 
the  PR  model.  The  natural  regression  strategy  (Fig.  7c) 
provides  the  poorest  resultant  threat  scores  and  these  reach 
their  peak  with  the  inclusion  of  the  fifth  predictor.  The 
predictors  selected  for  this  case  are  BMDl,  ENTR,  DVDP , 

PS  and  PBLD . 

The  functional  dependence,  A1/A0  statistics  and  96%/ 
05%  confidence  interval  values  for  this  model  can  be  seen  in 
Fig.  8.  As  with  the  PR  model,  the  specific  functional 
dependence  values  associated  with  the  selected  predictors 
never  fall  below  the  calculated  96%  functional  dependence 
confidence  interval. 

Figs.  9  and  10  show  the  contingency  tables  results 
for  the  EVAR  and  QUAD  threshold  models.  For  each  of  these 
models  the  independent  adjusted  VISCAT  I  threat  scores  have 
identical  values  of  0.32  (unadjusted,  0.38). 

The  two-stage  linear  regression  sequence  employed 
for  both  of  these  threshold  models  yields  very  similar  basic 
statistics.  For  the  EVAR  model,  a  threshold  value  of 
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0.648497  was  calculated  for  the  first-stage  VISCAT  I  versus 
VISCAT  II+III  separation.  This  threshold  was  based  upon  a 
VISCAT  I  sample  size  of  190  observations,  a  mean  of  0.659 
and  standard  deviation  of  0.205  and  a  combined  VISCAT  II+III 
sample  size  of  1722  observations,  a  mean  of  0.927  and 
standard  deviation  of  0.122.  The  second-stage  VISCAT  II 
versus  VISCAT  III  separation  was  based  on  a  calculated 
threshold  of  0.580128.  Associated  with  this  threshold  value 
were  311  VISCAT  II  observations  with  a  mean  of  0.708  and 
standard  deviation  of  0.142  and  1473  VISCAT  III  observations 
with  a  mean  of  0.850  and  standard  deviation  of  0.131. 

For  the  QUAD  model,  a  threshold  value  of  0.642104 
was  calculated  for  the  VISCAT  I  versus  VISCAT  II+III  first- 
stage  separation,  based  upon  the  sample  addressed  above.  A 
second-stage  quadratic  threshold  separating  the  VISCAT  II 
and  VISCAT  III  samples  was  calculated  to  be  0.580569.  This 
VISCAT  II  sample  contained  358  observations  with  a  mean  of 
0.643  and  standard  deviation  of  0.142  while  the  VISCAT  III 
sample  contained  1402  observations  with  a  mean  of  0.846  and 
standard  deviation  of  0.140. 

While  no  significant  difference  appears  to  exist 
between  the  results  of  the  two  threshold  models,  the  QUAD 
model  yields  a  slightly  higher  A0  and  slightly  lower  Al 
values  for  both  the  dependent  and  independent  data  sets. 

Table  V  shows  a  synopsis  of  the  key  statistical 
results  for  this  case.  The  best  models,  as  determined  by 


independent  adjusted  VISCAT  I  threat  scores,  are  the  two 
threshold  models.  Of  these  two  models,  the  QUAD  model 
achieves  the  highest  adjusted  AO,  namely  3.16%  (unadjusted, 
80.73%) . 

2.  Area  2,  TAU-24 

Fig.  11  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  AO  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR  model.  For  this  model,  adjusted  dependent  AO  values  of 
-0.03  and  adjusted  independent  threat  score  of  -.01  were 
obtained  for  grouping  sizes  four  through  nine.  At  the 
grouping  size  of  ten,  a  jump  in  scores  was  realized  and 
thus  ten  is  identified  as  the  only  possible  selection.  An 
associated  difficulty  in  utilizing  a  grouping  size  of  eight, 
nine  or  ten,  is  that  local  computer  storage  resources  are 
limited  to  two  megabytes.  This  decreases  the  usual  five 
predictor  array  to  only  four  predictors  as  witnessed  in  this 
case.  The  results  of  the  three  Preisendorf er  strategies  are 
shown  in  Figs.  12a  through  12c.  For  this  model,  the  MAXPROB 
I  and  MAXPROB  II  strategies  yield  identical  maximum  independent 
adjusted  VISCAT  I  threat  scores  of  0.21  (unadjusted,  0.27). 

For  each  of  the  maximum  probability  strategies,  an  initial 
threat  score  of  0.19  (unadjusted,  0.25)  was  achieved  with 
the  first  predictor,  E850,  solely.  The  slight  increase  to 
the  overall  peak  threat  score  was  obtained  with  the  inclusion 
of  the  second  predictor,  ENTR,  with  subsequent  independent 


VISCAT  I  threat  scores  decreasing  at  the  third  and  fourth 
predictor  levels.  Of  the  three  strategies,  natural  regression 
(Fig.  12c)  yielded  the  poorest  overall  threat  score  and  per¬ 
cent  correct  values.  These  relative  peak,  scores  for  the 
natural  regression  strategy  occur  with  the  inclusion  of  the 
fourth  and  final  predictor.  The  predictors  selected  for  this 
model  were:  E850,  ENTR,  DVDP  and  DIV925. 

The  associated  functional  dependence,  A0/A1  statis¬ 
tics  and  96%/05%  confidence  intervals  for  this  model  are 
shown  in  Fig.  13.  For  this  case,  the  third  and  fourth  pre¬ 
dictors'  root-sum-square  functional  dependence  values  exceed 
the  associated  96%  confidence  interval  values,  indicating 
significant  statistical  interdependence  of  these  predictors 
at  this  confidence  interval  level. 

Fig.  14  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  AO  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR+BMD  model.  The  dramatic  increase  in  independent  threat 
score  at  grouping  size  of  seven  identifies  it  as  the  optimal 
selection.  The  results  of  the  three  Preisendor fer  strategies 
are  shown  in  Figs.  15a  through  15c.  For  this  model,  the 
MAXPROB  I  and  MAXPROB  II  strategies  yield  identical  maximum 
independent  adjusted  VISCAT  I  threat  scores  of  0.26  (unad¬ 
justed,  0.32).  This  peak  score  was  achieved  with  the  inclu¬ 
sion  of  the  first  predictor.  In  this  case,  the  first  selected 
predictor  is  the  second  generated  linear  regression  equation 


44 


predictor,  BMD2  (Appendix  D) .  Following  the  initial  threat 
score  maxima,  the  scores  decreased  with  the  addition  of  the 
subsequent  four  predictors.  While  some  fluctuation  in  the 
threat  score  trend  was  observed  with  the  MAXPROB  II  strategy, 
independent  VISCAT  I  threat  scores  never  surpassed  their 
initial  maximum  value.  Of  the  three  strategies,  natural 
regression  (Fig.  15c)  provides  the  poorest  overall  indepen¬ 
dent  VISCAT  I  threat  score  of  0.21  (unadjusted,  0.28).  This 
score  was  achieved  with  the  addition  of  the  fifth  and  final 
predictor.  The  predictors  selected  for  this  model  were: 

BMD2,  VRT925 ,  ENTR,  U1000  and  RH. 

The  associated  functional  dependence,  A0/A1  statis¬ 
tics  and  96%/05%  confidence  intervals  are  shown  in  Fig.  16. 
For  this  model,  a  comparison  of  functional  dependence  and 
functional  dependence  96%  confidence  interval  values  indi¬ 
cates  that  the  final  three  predictors  have  root-sum-square 
functional  dependence  values  which  are  too  large  to  ensure 
significant  statistical  independence  at  the  96%  confidence 
interval  level. 

Figs.  17  and  18  show  the  contingency  tables  and 
associated  statistics  for  the  EVAR  and  QUAD  threshold  models. 
For  each  of  the  models,  the  independent  adjusted  VISCAT  I 
threat  scores  have  identical  values,  namely  0.29  (unadjusted, 
0.24).  The  two-stage  linear  regression  sequence  employed 
for  both  of  these  models  yields  fairly  similar  statistical 
results.  For  the  EVAR  model,  a  threshold  value  of  0.674932 


was  calculated  for  the  first-stage  VISCAT  I  versus  VISCAT 
II+III  separation  based  upon  a  VISCAT  I  sample  size  of  180, 
a  mean  of  0.682  and  a  standard  deviation  of  0.227  and  a 
VISCAT  II+III  sample  size  of  1580,  a  mean  of  0.938  and  a 
standard  deviation  of  0.109.  The  second-stage  VISCAT  II 
versus  VISCAT  III  separation  was  based  upon  a  calculated 
threshold  value  of  0.601717.  Associated  with  this  threshold 
were  300  VISCAT  II  observations  with  a  mean  of  0.733  and 
standard  deviation  of  0.149  and  1339  VISCAT  III  observations 
with  a  mean  of  0.857  and  standard  deviation  of  0.121. 

For  the  QUAD  model,  a  threshold  value  of  0.675210 
was  calculated  for  the  first-stage  VISCAT  I  versus  VISCAT 
II+III  separation  based  upon  the  sample  statistics  addressed 
above.  The  second-stage  threshold  separating  the  VISCAT  II 
and  VISCAT  III  samples  was  calculated  to  be  0.617455.  The 
VISCAT  II  sample  contained  300  observations  with  a  mean  of 
0.739  and  a  standard  deviation  of  0.125.  The  VISCAT  III 
sample  contained  1339  observations  with  a  mean  of  0.885  and 
standard  deviation  of  0.118. 

While  the  VISCAT  I  threat  scores  for  both  the  dependent 
and  independent  data  sets  are  identical  for  the  two  models, 
differences  in  other  statistics  are  apparent.  The  EVAR  model 
(Fig.  17),  for  example,  has  the  higher  independent  adjusted 
A0  scores,  namely  2.96%  (unadjusted,  81.34%),  as  compared  to 
scores  of  -63.31%  (unadjusted,  68.60%)  for  the  QUAD  model 
(Fig.  18)  . 
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In  general,  for  area  2,  TAU-24,  the  threshold  models 
again  provide  the  highest  independent  VISCAT  I  threat  scores 
(Table  V) .  Of  the  two  threshold  models,  the  EVAR  model  has 
a  slight  edge  in  AO  scores. 

3 .  Area  2,  TAU-48 

Fig.  19  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  AO  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR  model.  For  this  model,  the  initial  peak  values  of  dependent 
AO  and  independent  VISCAT  I  threat  score  at  the  grouping 
size  of  four  did  not  sufficiently  ascertain  four  as  the 
optimal  grouping  selection.  For  this  grouping  size,  the 
second  selected  predictor  ENTR  had  a  functional  dependence 
of  0.2952  as  compared  to  the  calculated  functional  dependence 
96%  confidence  interval  value  of  0.1932.  The  large  dis¬ 
parity  between  the  two  functional  dependence  values  indicates 
a  significant  statistical  correlation  between  E850  and  ENTR 
at  a  grouping  size  of  four  and  thus  grouping  size  four  was 
dropped  from  consideration.  The  selected  grouping  size  of 
nine,  which  unfortunately  carries  with  it  the  requirement 
of  a  very  large  computer  storage  forecast  array  at  the  fifth 
predictor  level,  had  a  functional  dependence  value  0.0930 
as  compared  to  a  functional  dependence  96%  confidence  interval 
value  of  0.0970  and  thus  was  selected  as  the  optimal  grouping 
size.  The  associated  functional  dependence,  A0/A1  statistics 
and  96%  confidence  intervals  are  shown  in  Fig.  20.  The  first 


three  predictors  selected  have  functional  dependence  values 
sufficiently  low  enough  to  ensure  no  significant  predictor 
interdependence . 

The  results  of  the  three  Preisendorf er  strategies 
are  shown  in  Figs.  21a  through  21c.  The  maximum  independent 
VISCAT  I  threat  score  achieved  for  the  three  strategies  was 
0.17  (unadjusted,  0.26)  and  was  obtained  with  the  MAXPROB  II 
strategy  with  the  addition  of  the  fifth  predictor.  It  should 
be  noted  that  the  independent  adjusted  VISCAT  I  threat 
scores  achieved  by  both  the  MAXPROB  I  and  MAXPROB  II  strate¬ 
gies  reached  near  peak  values  of  0.16  (unadjusted,  0.24) 
with  the  addition  of  the  second  predictor,  thus  greatly  mini¬ 
mizing  the  size  of  the  associated  forecast  array.  Of  the 
three  strategies,  natural  regression  (Fig.  21c)  yielded  the 
poorest  overall  adjusted  independent  VISCAT  I  threat  score, 
namely  0.09  (unadjusted,  0.18).  This  score  was  achieved 
with  the  inclusion  of  the  fourth  predictor  in  the  forecast 
array.  The  predictors  selected  for  this  model  were  E850, 
ENTR,  DVDP ,  DRAG  and  DIV925. 

Fig.  22  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  A0  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR+BMD  model.  For  this  model  a  grouping  size  of  nine  was 
selected.  The  results  of  the  MAXPROB  I,  MAXPROB  II  and 
natural  regression  strategies  are  shown  in  Figs.  23a  through 
23c.  For  this  model,  MAXPROB  I  and  MAXPROB  II  provide 


identical  maximum  independent  adjusted  VISCAT  I  threat  scores 
of  0.31  (unadjusted,  0.37).  These  scores  were  achieved  with 
the  inclusion  of  the  second  linear  regression  equation  pre¬ 
dictor  BMD2  (Appendix  D) .  For  each  of  these  strategies,  the 
independent  VISCAT  I  threat  scores  decrease  with  the  addi¬ 
tion  of  the  second  and  subsequent  predictors.  While  a  slight 
upward  progression  is  noticed  with  the  MAXPROB  II  strategy, 
the  peak  score  observed  at  the  first  predictor  level  is 
never  surpassed.  Of  the  three  Preisendorf er  strategies, 
natural  regression  (Fig.  23c),  yields  the  poorest  overall 
independent  VISCAT  I  threat  score,  namely  0.18  (unadjusted, 
0.26).  This  score  occurs  with  the  inclusion  of  the  fifth 
predictor  and  culminates  in  a  slow  increase  in  threat  score 
as  each  predictor  is  sequentially  added  to  the  forecast  array. 
The  predictors  selected  for  this  model  were  BMD2,  VRT925, 

ENTR,  U500  and  DRAG. 

Fig.  24  shows  the  functional  dependence,  A0/A1  sta¬ 
tistics  and  96%/05%  confidence  interval  values  for  the 
selected  predictors.  For  this  model,  the  second  and  third 
predictors'  functional  dependence  values  fall  below  the  96% 
confidence  interval  and  thus  are  not  significantly  inter¬ 
dependent  upon  one  another.  This  trend  changes  with  the 
fourth  and  fifth  predictors  which  have  functional  dependence 
values  greater  than  the  calculated  96%  confidence  interval 
values . 

Figs.  25  and  26  show  the  contingency  table  results 
for  the  EVAR  and  QUAD  threshold  models.  For  each  of  these 
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models,  the  independent  adjusted  VISCAT  I  threat  scores 
have  identical  values  of  0.21  (unadjusted,  0.29). 

The  two-stage  linear  regression  sequence  used  to 
separate  the  three  visibility  categories  yield  very  similar 
results  for  the  two  threshold  models.  For  the  EVAR  model, 
a  threshold  value  of  0.652554  was  calculated  for  the  first- 
stage  VISCAT  I  versus  VISCAT  II+III  sample  separation. 

This  threshold  value  is  based  upon  a  VISCAT  I  sample  size 
of  182  observations  with  a  mean  of  0.686  and  a  standard 
deviation  of  0.267  and  a  combined  VISCAT  II+III  sample  of 
1670  observations  with  an  associated  mean  of  0.930  and 
standard  deviation  of  0.106.  The  second  stage  VISCAT  II 
versus  VISCAT  III  regression  separation  yielded  a  threshold 
value  of  0.572257  based  upon  355  VISCAT  II  observations,  with 
a  mean  of  0.711  and  standard  deviation  0.135,  and  1408 
VISCAT  III  observations  with  a  mean  of  0.834  and  a  standard 
deviation  of  0.130. 

For  the  QUAD  model,  a  very  similar  threshold  value 
of  0.652554  was  calculated  for  the  first-stage  VISCAT  I 
versus  VISCAT  II+III  separation  based  upon  the  sample  first- 
stage  statistics  addressed  above.  A  second-stage  threshold 
value  of  0.564579  was  calculated  based  upon  330  VISCAT  II 
observations  with  a  mean  of  0.724  and  standard  deviation  of 
0.128,  and  1407  VISCAT  III  observations  with  a  mean  of  0.833 
and  a  standard  deviation  of  0.127. 


In  general,  the  results  of  these  two  threshold  models 


are  nearly  identical.  The  EVAR  model  shows  a  very  slight 
advantage  in  adjusted  independent  AO  scores,  namely  7.07% 
(unadusted,  80.11%)  as  compared  to  5.05%  (unadjusted, 

79.68%)  for  the  QUAD  model.  Similarly,  the  EVAR  model  yielded 
a  slightly  higher  independent  adjusted  threat  score  for  VISCAT 
I  combined  with  VISCAT  II  of  0.02  (unadjusted,  0.23)  versus 
an  adjusted  score  of  0.01  (unadjusted,  0.22)  for  the  QUAD 
model . 

For  Area  2,  TAU-48  the  PR+BMD  model  provides  the 
highest  overall  independent  VISCAT  I  threat  score  (Table  V) . 
The  difference  between  the  independent  adjusted  VISCAT  I 
threat  scores  for  the  PR+BMD  model  and  the  two  threshold 
models  is  minimal,  namely  0.02,  while  the  PR  model  is  0.14 
lower . 

B.  NORTH  ATLANTIC  OCEAN,  AREA  3W 

Area  3W  was  the  North  Atlantic  homogeneous  area  selected 
by  Karl  (1984)  for  his  initial  TAU-00  MOS  experimentation. 

This  area  borders  the  United  State's  eastern  seaboard  from 
the  vicinity  of  Cape  Charles,  Virginia  to  the  southeastern 
tip  of  Newfoundland.  The  area  encompasses  a  large  portion 
of  the  Georges  Banks  region  and  extends  to  approximately 
45°  W  longitude.  The  specific  detail  and  proximity  of  this 
area  can  be  seen  in  Fig.  2. 

Area  3W  constitutes  the  homogeneous  area  with  the  highest 
relative  frequency  of  VISCAT  I  observations  with  approximately 
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19%  of  the  total  number  of  visibility  observations  being  less 
than  2  kilometers  in  the  TAU-00,  TAU-24  and  TAU-48  periods. 

The  TAU-24  and  TAU-48  prognostic  periods  will  be  addressed 
in  this  document.  The  reader  is  advised  to  consult  Karl 
(1984)  for  detailed  information  concerning  area  3W,  TAU-00. 

1.  Area  3W,  TAU-24 

Fig.  27  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  AO  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR  model.  For  this  case  a  grouping  size  of  six  was  selected. 
Results  of  the  three  Preisendorf er  strategies  are  shown  in 
Figs.  28a  through  28c.  The  MAXPROB  II  strategy  achieves  a 
slightly  higher  independent  adjusted  VISCAT  I  threat  score 
of  0.21  (unadjusted,  0.36)  as  compared  to  a  score  of  0.20 
(unadjusted,  0.35)  for  the  MAXPROB  I  method.  For  each  of 
these  strategies,  the  maximum  threat  score  is  reached  with 
the  inclusion  of  the  fifth  and  final  predictor  in  the  fore¬ 
cast  array.  The  general  trend  of  these  two  strategies  is 
nearly  identical  and  show  an  initial  rise  in  threat  score 
at  the  first  predictor  level,  a  slight  decrease  with  the 
addition  of  the  second  and  third  predictors  and  a  secondary 
rise  at  the  fourth  and  fifth  predictor  levels.  The  poorest 
results  for  this  case  were  achieved  with  the  natural  regres¬ 
sion  strategy  (Fig.  28c) ,  for  which  an  independent  adjusted 
VISCAT  I  threat  score  of  0.16  (unadjusted,  0.32)  was  achieved. 
This  score  was  similarly  reached  with  the  addition  of  the 
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fifth  and  final  predictor.  The  predictors  selected  for  this 
model  were  DTDP ,  SHWRS,  ENTR,  U1000  and  DUDP . 

Fig.  29  shows  the  functional  dependence,  AO/Al 
statistics  and  96%/05%  confidence  intervals  for  this  model. 

For  this  case  only  the  second  predictor  has  a  functional 
dependence  value  which  falls  below  the  corresponding  96% 
confidence  interval  and  thus  meets  the  requisite  conditions 
regarding  predictor  interdependence.  Consequently,  the 
greatest  independent  threat  score  achieved,  which  coinci- 
dently  meets  the  functional  dependence  criteria,  occurs  with 
the  MAXPROB  II  strategy  at  the  inclusion  of  the  second  pre¬ 
dictor.  The  threat  score  achieved  in  this  particular  instance 
has  a  value  of  0.13  (unadjusted,  0.30). 

Fig.  30  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  AO  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR+BMD  model.  For  this  case  a  grouping  size  of  five  was 
selected.  Results  of  the  MAXPROB  I,  MAXPROB  II  and  natural 
regression  strategies  are  shown  in  Figs.  31a  through  31c. 

For  this  model,  the  two  maximum  probability  strategies  pro¬ 
vide  identical  peak  independent  adjusted  VISCAT  I  threat 
scores  of  0.28  (unadjusted,  0.42)  at  the  first  predictor 
level.  For  both  of  these  strategies,  the  addition  of  subse¬ 
quent  predictors  produces  a  steady  drop  off  in  threat  score 
values.  The  poorest  overall  results  for  this  case  are 
achieved  with  the  natural  regression  strategy  (Fig.  31c). 


This  method  yields  an  independent  adjusted  VISCAT  I  threat 
score  of  0.17  (unadjusted,  0.33)  which  was  obtained  with  the 
addition  of  the  fifth  and  final  predictor.  The  predictors 
selected  for  this  model  were  BMDl,  D500,  DVDP ,  ENTR  and 
U850  . 

Fig.  32  shows  the  associated  functional  dependence, 
A1/A0  statistics  and  96%/05%  confidence  interval  values  for 
the  predictors  chosen  for  this  model.  The  functional  depen¬ 
dence  versus  the  96%  confidence  interval  follows  a  peculiar 
trend  where  the  second  predictor  is  significantly  dependent 
upon  the  first  predictor  but  the  third  and  fourth  predictors 
are  conversely  sufficiently  uncorrelated  with  the  prior 
predictors  to  ensure  no  significant  functional  dependence. 

The  final  predictor  returns  to  being  functionally  dependent 
upon  the  previous  predictors.  This  trend  indicates  that  the 
relative  contribution  of  the  second  and  subsequent  predictors 
is  statistically  not  significant  at  the  preselected  96% 
confidence  interval  level. 

Figs.  33  and  34  show  the  contingency  table  results 
for  the  EVAR  and  QUAD  threshold  models.  The  results  of  these 
models  are  very  similar  with  the  EVAR  model  yielding  an 
independent  adjusted  VISCAT  I  threat  score  of  0.17  (unad¬ 
justed,  0.33)  as  compared  to  a  corresponding  threat  score 
of  0.16  (unadjusted,  0.32)  for  the  QUAD  model. 

For  the  EVAR  model,  a  first-stage  threshold  value  of 
0.561855  was  calculated  based  upon  270  VISCAT  I  observations 
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with  a  mean  of  0.590  and  standard  deviation  of  0.203  and 
1145  VISCAT  II+III  observations  with  a  mean  of  0.861  and 
standard  deviation  of  0.168.  The  second-stage  VISCAT  II 
versus  VISCAT  III  separation  was  based  upon  a  calculated 
threshold  value  of  0.542363.  Associated  with  this  threshold 
were  299  VISCAT  II  observations  with  a  mean  of  0.647  and 
standard  deviation  of  0.146  and  938  VISCAT  III  observations 
with  a  mean  of  0.794  and  standard  deviation  of  0.153. 

For  the  QUAD  model,  a  similar  threshold  value  of 
0.5559971  was  calculated  based  upon  the  first-stage  regres¬ 
sion  separation  listed  above.  A  second-stage  threshold 
value  of  0.540874,  separating  VISCAT  II  from  VISCAT  III,  was 
calculated  based  upon  305  VISCAT  II  observations  with  a  mean 
of  0.639  and  standard  deviation  of  0.157  and  940  VISCAT  III 
observations  with  a  mean  of  0.793  and  standard  deviation 
of  0.154. 

In  general,  the  PR+BMD  model  produced  the  best  overall 
results  for  this  case,  followed  by  the  PR  model  and  lastly 
the  two  threshold  models  (Table  V) .  The  independent  adjusted 
A0  score  of  the  PR+BMD  model,  which  corresponds  to  the  maxi¬ 
mum  independent  adjusted  VISCAT  I  threat  score,  is  similarly 
a  maximum  value  for  this  case,  namely  21.68%  (unadjusted, 
74.96%)  . 

2.  Area  3W,  TAU-48 

Fig.  35  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  A0  (dependent  data)  and  the 
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adjusted  VISCAT  I  threat  score  (independent  data)  for  the  PR 
model.  For  this  model,  an  equally  populous  grouping  size  of 
six  was  selected.  The  results  of  the  three  Preisendorf er 
strategies  are  shown  in  Figs.  36a  through  36c.  For  this 
case,  the  MAXPROB  II  strategy  achieves  the  highest  indepen¬ 
dent  adjusted  VISCAT  I  threat  score  of  0.18  (unadjusted, 

0.33)  as  compared  to  0.17  (unadjusted,  0.32)  for  the  MAXPROB 
I  strategy  and  0.12  (unadjusted,  0.22)  for  natural  regression. 
The  maximum  score  for  each  of  the  three  methods  was  achieved 
with  the  addition  of  the  fifth  and  final  predictor.  The 
statistical  score  trends  for  the  two  maximum  probability 
strategies  are  very  similar  and  reach  identical  near  peak 
independent  VISCAT  I  threat  scores  of  0.15  (unadjusted,  0.30) 
at  the  first  predictor  level.  This  is  particularly  note¬ 
worthy  when  considering  that  the  computer  forecast  array 
size  may  be  of  significant  operational  concern.  The  poorest 
strategy  for  this  case  is  natural  regression.  The  predictors 
selected  for  this  case  are  DTDP ,  SHWRS ,  ENTR,  U850  and 
DIV925. 

Fig.  37  shows  the  functional  dependence,  A0/A1  sta¬ 
tistics  and  96%/05%  confidence  interval  values  for  this 
model.  In  this  case,  only  the  second  predictor  strictly 
meets  the  requisite  functional  dependence  criteria  ensuring 
no  significant  dependence  of  one  predictor  upon  another.  The 
MAXPROB  II  independent  adjusted  A0  score,  which  corresponds 
to  the  peak  independent  VISCAT  I  threat  score  for  this  case, 


is  -2.47%  (unadjusted,  66.49%)  as  compared  to  AO  scores  of 
4.94%  (unadjusted,  68.91%)  for  MAXPROB  I  and  -16.87%  (unad¬ 
justed,  61.78%)  for  natural  regression. 

Fig.  38  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  AO  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR+BMD  model.  For  this  case  a  grouping  size  of  five  was 
selected.  The  results  of  the  three  Preisendorf er  strategies 
are  shown  in  Figs.  39a  through  39c.  For  this  case,  the 
MAXPROB  II  strategy  provides  the  highest  independent  ad¬ 
justed  VISCAT  I  threat  score  of  0.30  (unadjusted,  0.43). 

This  peak  score  slightly  surpasses  the  score  of  0.29 
(unadjusted,  0.42)  achieved  by  the  MAXPROB  I  method.  The 
trends  for  these  two  strategies  are  nearly  identical,  showing 
only  a  slight  oscillation  in  independent  threat  scores  as 
predictors  are  added.  The  peak  score  achieved  by  the  MAXPROB 
II  scheme  is  at  the  fifth  predictor  level  while  the  peak 
value  for  MAXPROB  I  is  obtained  with  the  inclusion  of  the 
first  predictor.  It  should  be  noted  that  the  results  at 
the  first  predictor  level  for  the  two  maximum  probability 
strategies  are  identical.  A  forecast  array  predicated  upon 
a  one  predictor  versus  five  predictor  array  size  requires 
four  orders  of  magnitude  less  computer  storage  resources  and 
is  therefore  a  desirable  characteristic  for  an  operational 
forecast  system.  Additionally,  the  independent  adjusted 
A0  scores,  achieved  by  both  schemes,  have  identical  maximum 


values  of  19.25%  (unadjusted,  73.76%)  at  the  first  predictor 
level  as  compared  to  a  maximum  value  of  12.76%  (unadjusted, 
71.47%)  for  the  natural  regression  strategy  at  the  fifth 
predictor  level.  The  poorest  strategy  for  this  case  is 
natural  regression  (Fig.  39c).  The  independent  VISCAT  I 
threat  scores  for  this  scheme  initially  yield  very  low  threat 
score  values  at  the  first  and  second  predictor  levels  with  a 
subsequent  rapid  rise  at  the  third,  fourth  and  fifth  predic¬ 
tor  levels.  This  rapid  rise  however  produces  a  threat  score 
value  of  only  0.19  (unadjusted,  0.34)  and  a  corresponding  AO 
of  -1.65%  (unadjusted,  66.76%)  at  the  fifth  and  final  predictor 
level.  The  predictors  selected  for  this  model  are  BMD2 , 

U1000,  ENTR,  DVDP  and  EAIR. 

Fig.  40  shows  the  functional  dependence,  A0/A1  sta¬ 
tistics  and  96%/0 5%  confidence  interval  values  for  this  model. 
For  this  case,  three  of  the  five  selected  predictors  do  not 
meet  the  96%  confidence  interval  criteria  for  functional 
independence.  This  further  justifies  the  use  of  a  single 
predictor  forecast  array  for  possible  operational  use. 

Figs.  41  and  42  show  the  contingency  table  results 
for  the  EVAR  and  QUAD  threshold  models.  The  results  of  these 
two  models  are  very  similar  with  the  EVAR  model  showing  a 
slight  advantage  in  independent  adjusted  VISCAT  I  threat 
score  of  0.15  (unadjusted,  0.33)  versus  0.14  (unadjusted, 

0.31)  for  the  QUAD  model.  Similarly,  the  EVAR  model  achieves 
a  slightly  higher  independent  adjusted  A0  of  13.17% 


(unadjusted,  71.60%)  versus  12.76%  (unadjusted,  71.47%) 
for  the  QUAD  model. 

For  the  EVAR  model,  a  first-stage  regression  thres¬ 
hold  value  of  0.577452  was  calculated  based  upon  a  VISCAT  I 
sample  size  of  290  observations  with  a  mean  of  0.620  and 
standard  deviation  of  0.211  and  a  combined  VISCAT  II+III 
sample  size  of  1197  observations  with  a  mean  of  0.860  and 
standard  deviation  of  0.153.  A  second-stage  threshold  value 
of  0.548587  separating  VISCAT  II  and  VISCAT  III  was  calcu¬ 
lated  based  upon  a  VISCAT  II  sample  size  of  328  with  a  mean 
of  0.654  and  standard  deviation  of  0.142  and  971  VISCAT  III 
observations  with  a  mean  of  0.777  and  standard  deviation  of 
0.136. 

The  first-stage  threshold  value  of  0.572592  for  the 
QUAD  model  was  generated  with  the  above  VISCAT  I  versus 
VISCAT  II+III  sample  statistics.  A  second-stage  threshold 
value  of  0.548717  was  based  upon  333  VISCAT  II  observations 
with  a  mean  of  0.649  and  standard  deviation  of  0.138. 

In  general,  the  model  which  produces  the  highest 
independent  VISCAT  I  threat  score  for  this  case  is  the 
PR+BMD  model  while  the  highest  independent  A0  score  is 
achieved  with  the  EVAR  threshold  model  (Table  V) .  The  rela¬ 
tively  large  independent  threat  score  dominates  the  scores 
however,  and  therefore  the  PR+BMD  model  is  determined  to  be 
the  optimal  model  in  this  case. 


C.  NORTH  ATLANTIC  OCEAN,  AREA  4 

Area  4  was  selected  for  evaluation  because  of  its  rela¬ 
tively  low  frequency  (approximately  3%  of  the  total)  of 
VISCAT  I  observations.  It  was  hoped  that  this  area  would 
statistically  represent  a  region  where  there  was  an  insuffi¬ 
cient  number  of  VISCAT  I  observations  to  allow  for  study  of 
a  forecast  region  where  results  were  anticipated  to  be  poor, 
yet  enough  VISCAT  I  observations  to  avoid  any  "rare  event" 
statistical  entanglements. 

This  area  encompasses  a  broad  region  of  the  North  Atlantic 
Ocean  which  is  generally  to  the  south  of  area  2  and  east  and 
southeast  of  area  3W.  Area  4's  southern  border  reaches  to 
the  northeastern  tip  of  Portugal  and  extends  northward  through 
the  English  Channel  to  encompass  the  southern  portion  of  the 
North  Sea. 

1 .  Area  4,  TAU-00 

Fig.  43  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  AO  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR  model.  Several  unique  characteristics  were  encountered 
for  this  case  which  had  not  been  previously  been  observed. 

The  previously  observed  variation  of  dependent  AO  and  inde¬ 
pendent  threat  scores,  associated  with  the  sequential  varia¬ 
tion  in  grouping  size  from  four  through  ten,  was  not  initially 
achieved.  For  this  case,  non-zero  values  of  dependent  AO 
and  independent  VISCAT  I  threat  score  were  only  achieved 


after  three  iterations  of  the  predictor  selection  procedure. 

The  grouping  size  of  four  was  deleted  from  consideration 
in  the  third  iteration  because  the  associated  AO  value, 
achieved  at  that  predictor  level,  did  not  exceed  the  previous 
AO  value  at  the  second  predictor  level.  For  this  case,  the 
independent  VISCAT  I  threat  scores  maintained  indentically 
low  values,  while  a  relative  peak  in  AO  was  achieved  at  a 
grouping  size  of  eight.  For  this  reason,  eight  was  selected 
as  the  optimal  grouping  size  for  this  model. 

Figs.  44a  through  44c  represent  the  results  of  the 
three  Preisendorf er  strategies.  For  each  of  the  schemes, 
the  independent  VISCAT  I  threat  scores  at  the  first  three 
predictor  levels  reveal  the  near-zero  scores  encountered  in 
the  grouping  size  selection  process.  The  highest  independent 
adjusted  VISCAT  I  threat  score,  namely  0.08  (unadjusted, 

0.11)  is  achieved  with  the  MAXPROB  II  strategy  at  the  fifth 
and  final  predictor  level.  For  this  model,  the  MAXPROB  I 
and  natural  regression  strategies  yield  only  slightly  inferior, 
identical  independent  adjusted  threat  scores  of  0.04  (inad- 
justed,  0.07)  which  are  achieved  at  the  fifth  predictor 
level.  The  MAXPROB  I  strategy  yields  the  highest  independent 
adjusted  A0  score  of  -15.77%  (unadjusted,  82.45%)  as  compared 
to  scores  of  -28.63%  (unadjusted,  80.50%)  for  natural 
regression  and  -34.85%  (unadjusted,  79.56%)  for  the  MAXPROB 
II  strategy.  The  predictors  selected  for  this  model  are 
V500,  DVDP ,  STRTTH ,  E500  and  ENTR. 
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Fig.  45  shows  the  functional  dependence,  A0/A1  sta¬ 
tistics  and  96%/05%  confidence  interval  values  for  this  model. 
In  this  particular  case,  only  the  third  predictor  displays 
a  functional  dependence  value  less  than  the  96%  confidence 
interval  value.  This  renders  the  threat  scores  achieved  by 
this  model,  beyond  the  first  predictor  level,  statistically 
not  significant,  if  strict  adherence  to  the  basic  functional 
dependence  criteria  is  followed. 

Fig.  46  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  AO  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR+BMD  model.  For  this  case  a  grouping  size  of  nine  was 
selected.  The  results  of  the  three  Preisendorfer  strategies 
are  shown  in  Figs.  47a  through  47c.  Generally,  the  results 
for  this  model  differ  very  little  from  the  previously  dis¬ 
cussed  PR  model.  This  case  reflects  the  first  and  only 
occurrence  where  the  Preisendorfer  methodology  coupled  with 
linear  regression  equation  predictors  (PR+BMD  model)  did  not 
yield  superior  results  to  the  PR  model.  The  trends  for 
these  three  strategies  are  generally  quite  similar.  The 
MAXPROB  II  scheme  provides  the  highest  independent  adjusted 
VISCAT  I  threat  score  of  0.09  (unadjusted,  0.11),  as  com¬ 
pared  to  scores  of  0.08  (unadjusted,  0.11)  for  MAXPROB  I 
and  0.07  (unadjusted,  0.10)  for  natural  regression.  For  each 
of  these  three  strategies,  the  maximum  independent  VISCAT  I 
threat  score  was  achieved  with  the  inclusion  of  the  fifth 


and  final  predictor.  The  independent  AO  scores  associated 
with  the  peak  threat  scores  are  near  their  lowest  values  at 
the  fifth  predictor  level  with  the  MAXPROB  I  scheme  yielding 
the  highest  relative  independent  adjusted  AO  of  -19.09% 
(unadjusted,  81.95%)  followed  by  natural  regression  with  a 
score  of  -26.56%  (unadjusted,  80.82%)  and  lastly  MAXPROB  II 
with  a  score  of  -39.42%  (unadjusted,  78.87%).  The  predictors 
selected  for  this  model  are  BMD2 ,  DUDP,  ENTR ,  DEDP  and 
U1000. 

Fig.  48  shows  the  functional  dependence,  A0/A1  sta¬ 
tistics  and  96%/05%  confidence  interval  values  for  this 
model.  Generally,  the  relative  difference  between  the  func¬ 
tional  dependence  and  96%  functional  dependence  confidence 
interval  values  is  much  less  severe  than  with  the  previously 
discussed  model.  While  only  the  third  predictor's  functional 
dependence  value  meets  the  96%  confidence  interval  criteria 
for  significance,  the  other  predictors  are  only  marginally 
insignificant. 

The  application  of  the  EVAR  and  QUAD  threshold  models 
to  this  case  presented  results  which  had  not  been  previously 
encountered.  The  first-stage  VISCAT  I  versus  VISCAT  II+III 
separation  calculation  results  in  a  QUAD  threshold  value 
which  is  imaginary  and  an  unrealistic  EVAR  threshold  value 
of  209.588882.  These  thresholds  were  calculated  based  upon 
a  VISCAT  I  sample  size  of  85  observations  with  a  mean  of 
-1.012  and  standard  deviation  of  6.280  and  a  combined  VISCAT 
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II+III  sample  size  of  3096  observations  with  a  mean  of 
-1.864  and  standard  deviation  of  7.092.  These  results  are 
linked  to  the  preponderance  of  VISCAT  III  observations  in 
this  area  coupled  with  the  fact  that  these  employed  threshold 
models  are  designed  to  provide  for  a  minimum  error  when 
separating  samples.  These  results  indicate  that  a  forecast 
model  predicated  upon  the  dependent  data  set  employed  in  this 
case  would  strictly  forecast  VISCAT  III. 

2.  Area  4,  TAU-24 

Fig.  49  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  AO  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR  model.  This  case  required  three  iterations  of  the  four 
through  ten  grouping  size  calculations  before  any  non-zero 
dependent  AO  or  independent  VISCAT  I  threat  score  values  were 
achieved.  Additionally,  for  the  grouping  size  of  four,  no 
increase  in  AO  was  observed  at  the  second  predictor  level 
and  therefore  was  deleted  from  consideration.  A  grouping 
size  of  five  was  ultimately  selected  for  this  model. 

Figs.  50a  through  50c  represent  the  results  of  the 
three  Preisendorf er  strategies.  Generally,  the  independent 
VISCAT  I  threat  scores  yielded  for  these  schemes  are  poor 
with  the  highest  independent  adjusted  VISCAT  I  threat  score 
of  0.05  (unadjusted,  0.07)  being  achieved  by  the  MAXPROB  II 
strategy  at  the  fifth  predictor  level  followed  by  MAXPROB  I 
with  a  score  of  0.02  (unadjusted,  0.05)  and  natural  regression 


with  0.01  (unadjusted,  0.04).  The  AO  scores  corresponding 
to  these  values  provide  for  a  slightly  different  scoring 
hierarchy.  The  highest  independent  adjusted  AO  score, 
namely  0.94%  (unadjusted,  85.64%),  is  attained  by  the 
MAXPROB  I  strategy  as  compared  to  scores  of  -30.05%  (unad¬ 
justed,  81.34%)  for  natural  regression  and  -37.56%  (unad¬ 
justed,  80.05%)  for  MAXPROB  II.  The  predictors  selected  for 
this  model  are  VRT925,  DTDP ,  ENTR ,  V850  and  DVRTDP . 

Fig.  51  shows  the  functional  dependence,  A0/A1  sta¬ 
tistics  and  96%/05%  confidence  interval  values  for  this 
model.  For  this  case,  each  predictor  following  VRT925  proved 
to  be  significantly  functionally  dependent  on  its  predecessors 
and  therefore  only  a  single  predictor  forecast  array  is 
justifiable  for  this  model. 

Fig.  52  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  A0  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR+BMD  model.  As  in  the  previous  case,  three  iterations  of 
dependent  A0  and  independent  VISCAT  I  calculations  were 
required  before  any  non-zero  scores  were  achieved.  Addi¬ 
tionally,  in  this  case,  the  grouping  sizes  of  four  and  five 
were  deleted  from  consideration  as  they  did  not  provide  an 
increase  of  A0  at  the  second  predictor  levels.  The  grouping 
size  ultimately  selected  for  this  model  was  nine. 

Figs.  53a  through  53c  show  the  results  of  the  three 
Preisendorfer  strategies  for  this  model.  The  scores  for 
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this  model,  as  in  the  previously  described  case,  are  quite 
poor  and  show  very  little  improvement  over  the  PR  model. 

The  highest  independent  adjusted  VISCAT  I  threat  score, 
namely  0.06  (unadjusted,  0.09),  was  achieved  by  the  MAXPROB 
II  strategy  followed  by  scores  of  0.05  (unadjusted,  0.07) 
for  the  MAXPROB  I  strategy  and  0.03  (unadjusted,  0.06)  for 
natural  regression.  The  corresponding  independent  adjusted 
A0  scores  show  a  maximum  score  of  -19.25%  (unadjusted, 

82.71%)  for  the  MAXPROB  I  strategy  followed  by  scores  of 
-30.05%  (unadjusted,  81.14%)  for  natural  regression  and 
-39.91%  (unadjusted,  79.71%)  for  the  MAXPROB  II  strategy. 

Fig.  54  shows  the  functional  dependence,  A0/A1  sta¬ 
tistics  and  96%/05%  confidence  interval  values  for  this 
model.  For  this  case,  the  relative  magnitude  of  the  differ¬ 
ence  between  the  actual  functional  dependence  and  its  96% 
confidence  interval  value  is  quite  small.  It  is  only  at  the 
second  predictor  level,  that  the  calculated  values  do  not 
exceed  the  corresponding  96%  confidence  interval  value. 

Fig.  55  shows  the  contingency  table  results  for  the 
EVAR  threshold  model.  The  QUAD  model  provided  an  imaginary 
threshold  value  at  the  second  regression  stage  and  therefore 
did  not  allow  completion  of  the  entire  separation  sequence. 

This  represents  the  only  occurrence  where  a  valid  equal 
variance  threshold  was  calculated  but  a  corresponding  quadratic 
threshold  proved  to  be  imaginary. 
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The  results  of  the  EVAR  model  were  in  keeping  with 
those  of  the  previously  described  PR  and  PR+BMD  models.  An 
independent  adjusted  VISCAT  I  threat  score  of  0.05  (unad¬ 
justed,  0.07),  was  achieved  with  a  corresponding  independent 
adjusted  AO  value  of  -13.15%  (unadjusted,  83.59%). 

The  first-stage  regression  separation  for  this  model 
was  based  upon  a  calculated  threshold  value  of  0.908275. 
Associated  with  this  threshold  were  449  VISCAT  I  observations 
with  a  mean  of  0.953  and  standard  deviation  of  0.030  and 
2489  VISCAT  II+III  observations  with  a  mean  of  0.976  and 
standard  deviation  of  0.027.  The  second-stage  VISCAT  II 
versus  VISCAT  III  separation  was  based  upon  a  calculated 
threshold  value  of  0.683569.  Associated  with  this  threshold 
is  a  VISCAT  II  sample  size  of  69  observations  with  a  mean  of 
0.831  and  standard  deviation  of  0.066  and  887  VISCAT  III 
observations  with  a  mean  of  0.912  and  standard  deviation  of 
0.078. 

For  the  QUAD  model,  an  initial  first-stage  threshold 
value  of  0.908275  was  successfully  calculated  with  the 
sample  statistics  addressed  above.  The  second-stage  regres¬ 
sion  attempt  was  based  upon  a  VISCAT  II  sample  size  of  65 
observations  with  a  mean  of  0.829  and  standard  deviation  of 
0.067  and  VISCAT  III  sample  size  of  853  observations  with  a 
mean  of  0.905  and  standard  deviation  of  0.079.  These  sample 
statistics  produced  an  imaginary  threshold  value. 

In  general  Area  4,  TAU-24  is  characterized  by  very 
poor  independent  VISCAT  I  threat  scores.  This  indicates 
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that  there  is  very  little  skill  in  forecasting  visibility 
conditions  of  less  than  or  equal  to  2  kilometers  in  this 
area.  The  evaluated  models  show  little  variation  in  scores 
with  the  best  relative  model  for  this  area  and  forecast 
projection  being  the  PR+BMD  model  (Table  V) . 

3 .  Area  4,  TAU-48 

Fig.  56  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  AO  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the  PR 
model.  As  in  the  TAU-00  and  TAU-24  forecast  projections  for 
this  area,  the  calculation  and  evaluation  of  dependent  AO 
and  independent  VISCAT  I  threat  scores  had  to  be  run  through 
three  iterations  before  any  non-zero  statistics  were  obtained. 
The  grouping  sizes  of  four  and  five  were  deleted  from  con¬ 
sideration  because  the  addition  of  predictors  at  those  grouping 
sizes  did  not  provide  for  any  increase  in  AO  scores.  Based 
on  an  evaluation  of  the  results  as  shown  on  Fig.  56,  a  group¬ 
ing  size  of  seven  was  selected. 

Figs.  57a  through  57c  show  the  results  of  the  three 
Preisendorf er  strategies  for  this  model.  In  general,  the 
near-zero  statistical  scores  encountered  in  the  grouping 
selection  process,  can  be  seen  through  the  third  predictor 
level,  along  with  a  noticeable  increase  in  scores  at  the 
fourth  and  fifth  predictor  level.  The  MAXPROB  I  strategy 
yields  the  highest  independent  adjusted  VISCAT  I  threat 
score,  namely  0.18  (unadjusted,  0.20),  for  this  model 


followed  by  a  natural  regression  score  of  0.16  (unadjusted, 
0.19)  and  lastly  by  MAXPROB  II  with  a  score  of  0.13  (unad¬ 
justed,  0.16).  This  is  the  first  and  only  encountered  case 
where  the  natural  regression  strategy  effectively  achieved  a 
maximum  independent  VISCAT  I  threat  score  which  is  higher 
than  either  of  the  two  maximum  probability  strategies.  The 
independent  AO  scores  associated  with  these  peak  independent 
VISCAT  I  threat  scores,  adhere  to  this  same  scoring  sequence, 
with  MAXPROB  I  achieving  a  value  of  -14.94%  (unadjusted, 
81.86%)  followed  by  natural  regression  with  a  score  of 
-27.80%  (unadjusted,  79.83%)  and  MAXPROB  II  with  an  indepen¬ 
dent  adjusted  A0  of  -43.98%  (unadjusted,  77.78%). 

Fig.  58  shows  the  functional  dependence,  A0/A1  sta¬ 
tistics  and  96%/05%  confidence  interval  values  for  this  model. 
In  this  case,  only  the  third  predictor's  functional  dependence 
value  falls  below  the  associated  96%  confidence  interval 
value.  The  predictors  selected  for  this  model  are  VRT925, 
DVRTDP ,  ENTR ,  DUDP  and  RH. 

Fig.  59  shows  the  relationship  of  equally  populous 
grouping  size  to  the  adjusted  A0  (dependent  data)  and  the 
adjusted  VISCAT  I  threat  score  (independent  data)  for  the 
PR+BMD  model.  As  in  the  previous  area  4  cases,  three  com¬ 
plete  iterations  of  the  four  through  ten  grouping  size  calcu¬ 
lations  had  to  be  performed  before  any  non-zero  dependent 
A0  or  independent  VISCAT  I  values  were  achieved.  The 
grouping  size  ultimately  selected  for  this  model  was  nine. 


Figs.  60a  through  60c  represent  the  results  of  the 
three  Preisendor fer  strategies.  A  unique  result  of  this 
model  is  that  for  the  first  time,  the  PR+BMD  model  did  not 
achieve  independent  VISCAT  I  threat  scores  which  exceed  those 
achieved  by  the  PR  model.  The  peak  independent  adjusted 
VISCAT  I  threat  score  of  0.17  (unadjusted,  0.20)  is  achieved 
by  the  MAXPROB  II  strategy  at  the  third  predictor  level. 

The  predictors  selected  for  this  model  are  BMD1,  DDVDP, 

DUDP ,  ENTR  and  PRECIP. 

Fig.  61  shows  the  functional  dependence,  A0/A1  sta¬ 
tistics  and  96%/05%  confidence  interval  values  for  this 
model.  Only  the  second  predictor  sufficiently  meets  the  96% 
confidence  interval  significance  criteria.  Based  upon  a 
strict  adherence  to  the  preselected  96%  confidence  interval 
significance  requirements,  these  functional  dependence 
values  provide  cause  for  uncertainty  in  the  representative¬ 
ness  of  the  scores  achieved  after  the  second  predictor  level. 

Figs.  62  and  63  show  the  contingency  table  results 
for  the  EVAR  and  QUAD  threshold  models.  The  results  of  these 
models  are  very  similar  and  generally  quite  poor.  The  QUAD 
model  achieves  the  highest  relative  independent  adjusted 
VISCAT  I  threat  score  of  0.01  (unadjusted,  0.04)  as  compared 
to  a  score  of  -0.01  (unadjusted,  0.02)  for  the  EVAR  model. 

The  QUAD  model  similarly  achieves  the  highest  independent 
adjusted  A0  score  of  -2.07%  (unadjusted,  83.89%)  versus  a 
score  of  -7.88%  (unadjusted,  82.97%)  for  the  EVAR  model. 


For  the  EVAR  model,  a  first-stage  regression  thres¬ 
hold  value  of  0.847203  was  calculated  based  upon  a  VISCAT  I 
sample  size  of  109  observations  with  a  mean  of  0.918  and 
standard  deviation  of  0.052  and  2947  VISCAT  II+III  observa¬ 
tions  with  a  mean  of  0.967  and  standard  deviation  of  0.037. 

The  second-stage  VISCAT  II  versus  VISCAT  III  separation  was 
based  upon  a  calculated  threshold  value  of  0.629338.  Asso¬ 
ciated  with  this  threshold  is  495  VISCAT  II  observations  with 
mean  of  0.861  and  standard  deviation  of  0.103. 

For  the  QUAD  model,  a  first-stage  separation  thres¬ 
hold  value  of  0.847203  was  calculated  upon  the  associated 
sample  statistics  addressed  above.  A  second-stage  threshold 
value  of  0.613739  was  calculated  based  upon  481  VISCAT  II 
observations  with  a  mean  of  0.770  and  standard  deviation  of 
0.089  and  2522  VISCAT  III  observations  with  a  mean  of  0.862 
and  standard  deviation  of  0.100. 

The  overall  results  associated  with  the  area  4,  TAU- 
48  case  are  particularly  unique.  The  independent  adjusted 
TAU-48  VISCAT  I  threat  score  represents  the  highest  area  4 
independent  VISCAT  I  threat  score  (by  a  minimum  of  0.09) 
achieved,  as  compared  to  TAU-00  or  TAU-24.  The  maximum 
independent  VISCAT  I  threat  score  is  achieved  by  the  PR 
model . 

Following  the  completion  of  the  testing  and  evaluation 
of  the  FATJUNE  1983  data  set,  a  series  of  preliminary  experi¬ 
ments  were  performed  with  the  May  15  to  June  23  1984  data 
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set  for  the  TAU-24  model  forecast  projection.  These  experi¬ 
ments  consisted  of  evaluating  the  TAU-24,  1983  forecast 
arrays  and  equations  (generated  with  FATJUNE  1983  data) 
with  training  and  testing  cases  of  TAU-24,  1984  data.  In 
performing  this  evaluation,  the  1984  data  set  was  divided 
into  "dependent"  and  "independent"  portions.  This  data 
separation  is  a  function  of  the  specific  mechanics  of  the 
computer  programs  utilized  in  this  study  and  is  not  associated 
with  the  generation  of  additional  forecast  arrays  or  equations. 
Two  homogeneous  areas  were  evaluated,  namely  area  2  and  area 
3W.  This  essentially  provided  an  independent  verification 
of  the  utility  of  the  1983  forecast  arrays  and  equations  in 
predicting  observed  1984  visibility  in  these  areas. 

In  general,  the  skill  and  contingency  table  results 
for  these  experiments  compare  very  favorably  to  those  achieved 
with  the  FATJUNE  1983  data.  A  summary  of  the  results  of 
each  of  the  evaluated  models  is  provided  in  Table  VII.  For 
area  2,  a  peak  independent  adjusted  VISCAT  I  threat  score 
(1984  data),  namely  0.27  (0.33  unadjusted),  was  achieved 
with  each  of  the  two  threshold  models.  This  compares  to  a 
peak  independent  adjusted  VISCAT  I  threat  score  (1983  data) 
of  0.29  (0.34  unadjusted)  achieved  by  the  same  two  models. 

For  area  3W,  a  peak  independent  adjusted  threat  score  (1984 
data)  of  0.28  (unadjusted,  0.42)  similarly  compares  to  a 
peak  independent  adjusted  threat  score  (1983)  of  0.13 
(unadjusted,  0.36). 


V.-  ^  ^  ^  ^  ^  >.1  ,■»  O  C.  V  V  V.  U  V.  V.  V.  V 
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The  overall  results  of  the  1984  data  experiments  can 
be  seen  in  Table  VII  and  are  represented  by  Fig.  64  which 
illustrates  the  results  of  the  PR+BMD  model,  MAXPROB  II 
strategy,  for  area  2,  TAU-24. 

A  review  of  the  results  associated  with  area  4  for 
TAU-00,  TAU-24  and  TAU-48  indicates  that  none  of  the  models 
evaluated  achieved  very  encouraging  skill  and  threat  scores. 
Consequently,  the  maximum-likelihood-of -detection  criteria 
(MLDC)  was  proposed  as  an  alternative  technique  to  increase 
threat  scores  in  area  4. 

A  series  of  experiments  involving  an  arbitrarily 
selected  population  of  two  hundred  normally  distributed 
events,  partitioned  into  eight  separate  threat/non- threat 
samples,  were  performed  to  demonstrate  the  theoretical  utility 
of  the  MLDC  at  low  threat  frequencies.  Threshold  values 
were  calculated,  for  various  threat  frequencies,  using  the 
EVAR  minimum  probable  error  and  MLDC  techniques  and  two  by 
two  contingency  tables  were  constructed  to  tabulate  the  asso¬ 
ciated  threat  score,  percent  correct  and  false  alarm  rate 
results.  Fig.  65  shows  the  resulting  plot  of  threat  score 
versus  threat  frequency  which  illustrates  the  amount  of  in¬ 
crease  in  threat  score  associated  with  the  MLDC  model.  Asso¬ 
ciated  with  these  higher  threat  scores  are  correspondingly 
higher  "costs,"  namely  higher  false  alarm  rates,  illustrated 
in  Fig.  66,  and  lower  percent  correct  scores,  illustrated  in 


A  set  of  two  experiments  was  performed,  utilizing 
FATJUNE  1983,  TAU-24  data  and  the  two-stage  separation 
scheme  outlined  in  Chapter  IV  (MLDC  model) ,  to  evaluate  the 
relative  performance  of  the  MLDC  and  EVAR  models  on  area  4 . 

In  general,  the  results  of  these  two  experiments  were  consis¬ 
tent  with  the  results  predicted  by  the  aforementioned 
theoretical  experiments.  The  most  obvious  area  of  agreement 
is  the  significantly  lower  independent  adjusted  VISCAT  I 
and  VISCAT  II  threat  scores  (both  are  considered  threatening 
events  in  this  study),  namely  0.01  (unadjusted,  0.04)  and 
-0.14  (unadjusted,  0.00),  achieved  by  the  EVAR  model.  Fig. 

68,  as  compared  to  the  corresponding  scores  of  0.04  (unad¬ 
justed,  0.07)  and  0.03  (unadjusted,  0.15)  achieved  by  the 
MLDC  model.  Fig.  69. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  primary  objective  of  this  study  was  to  expand  upon 
the  initial  research  and  experimentation  presented  by  Karl 
(1984)  and  to  propose  a  viable  statistical  forecasting  scheme 
suitable  for  eventual  employment  in  an  operational  U.S.  Navy 
marine  visibility  MOS  forecasting  system.  In  general, 
while  the  results  of  linear  regression  and  the  evaluated 
Preisendorf er  models  are  roughly  comparable,  it  has  been 
shown  that  two  specific  statistical  approaches,  namely  the 
PR+BMD  model's  MAXPROB  II  strategy  and  the  linear  regression 
models,  yield  the  best  results  (as  measured  by  independent 
VISCAT  I  threat  score)  achieved  in  this  study.  The  PR+BMD 
model  achieved  the  best  results  for  six  of  the  eight  evaluated 
cases:  area  2,  TAU-48;  area  3W,  TAU-24  and  TAU-48;  and  area 

4,  TAU-00 ,  TAU-24  and  TAU-48.  The  nearly  identical  results 
of  both  the  equal  variance  and  quadratic  linear  regression 
threshold  models  provided  the  best  skill  and  threat  scores 
for  area  2,  TAU-00  and  TAU-24.  A  common  characteristic  of 
each  of  the  evaluated  cases  is  that  the  predictability  of 
visibility  category  II  is  relatively  very  poor  and  nearly 
always  poorer  than  that  for  visibility  categories  I  or  III. 
This  pattern  affirms  the  findings  of  similar  Pacific  Ocean 
visibility  studies  [Renard  and  Thompson,  1984]  as  well  as 
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those  documented  by  Karl  (1984)  and  further  supports  Karl's 
recommendation  to  change  from  a  three-category  to  a  two- 
category  visibility  forecasting  scheme. 

An  evaluation  of  the  overall  results  of  this  study  shows 
that  no  real  connection  between  individual  model/strategy 
and  either  the  homogeneous  oceanic  area  (2,  3W  and  4)  or 
model  output  time  (TAU-00,  TAU-24  and  TAU-48)  can  be  made. 

The  linear  regression  threshold  models  performed  best  for 
area  2,  the  intermediate  poor  visibility  oceanic  area,  while 
the  Preisendorfer  methodology  incorporating  linear  regression 
equation  predictors  proved  the  best  in  the  evaluated  homogene¬ 
ous  areas  with  the  greatest  and  lowest  relative  concentration 
of  poor  visibility  observations.  The  trend  of  visibility 
category  I  skill  and  threat  scores,  for  each  homogeneous  area 
and  model  output  time,  seems  to  contradict  the  preliminary 
supposition  that  peak  skill  scores  would  be  associated  with 
the  area  containing  the  greatest  frequency  of  poor  visiblity 
observations  and  the  TAU-00  model  output  time.  This  result 
is  most  apparent  with  area  4,  where  threat  scores  increase 
with  increasing  model  forecast  projections  until  they  achieve 
values  at  TAU-48,  which  are  nearly  identical  to  those  for  the 
other  two  homogeneous  areas.  This  type  of  trend  in  skill 
and  threat  scores  most  likely  reflects  the  overall  strength 
of  the  statistical  relationships  for  the  predictand/predictors 
involved  irrespective  of  the  frequency  of  specific  visibility 


observations . 


In  several  cases,  the  maximum  independent  visiblity 
category  I  threat  score  achieved  by  the  PR+BMD  model  was 
reached  at  the  first  predictor  level.  In  several  additional 
cases,  threat  score  values  which  were  only  marginally  lower 
than  peak  value,  were  similarly  achieved  at  the  first  pre¬ 
dictor  level.  Forecasting  arrays  involving  only  one  predic¬ 
tor  drastically  reduce  required  computer  storage  and 
consequently  such  arrays  are  a  desirable  attribute  to  any 
operational  MOS  forecasting  system.  A  MOS-type  forecasting 
system  predicated upon  such  a  small  number  of  predictors 
would  prove  extremely  beneficial  in  an  independent  single 
station  forecasting  scenario  such  as  that  experienced  by  an 
aircraft  carrier  based  U.S.  Navy  Oceanography  Officer. 

The  concept  and  practical  employment  of  functional 
dependence,  associated  with  the  Preisendorf er  methodology, 
provides  a  greater  restriction  on  the  statistical  significance 
of  the  skill  and  threat  score  results  achieved  in  this  study, 
as  compared  to  that  which  was  previously  experienced  by 
Karl  (1984) .  it  was  shown  that  the  calculated  functional 
dependence  values  for  each  respective  predictor  or  group  of 
predictors  often  exceeded  the  associated  96%  confidence 
interval  value  at  the  first  or  second  predictor  level  and 
rarely  met  the  requirements  for  significance  for  an  entire 
array  of  selected  predictors.  This  restriction  further 
indicates  that  any  operational  forecasting  scheme  should 
most  likely  be  composed  of  only  a  minimal  number  of  select 
predictors . 


The  difference  between  the  equal  variance  and  quadratic 
threshold  models  was  shown  to  be  very  minimal.  The  two-stage 
visibility  category  separation  approach  is  designed  to 
handle  cases  with  distinct  separability  between  categories 
while  providing  for  minimal  error  in  the  calculated  threshold 
values.  This  condition  was  not  met  in  the  area  4,  TAU-00 
and  TAU-24  cases  and  subsequently  lead  to  unrealistic  thres¬ 
hold  values. 

The  preliminary  independent  evaluation  of  the  15  May  to 
23  June  1984  data  set  provided  a  crucial  test  and  verifica¬ 
tion  of  the  utility  of  the  forecast  arrays  and  equations 
presented  in  this  study. 

The  introduction  and  initial  evaluation  of  the  maximum- 
likelihood-of-detection  threshold  model  offers  another 
technique  to  the  pool  of  visibility  prediction  schemes. 

This  method  appears  to  be  most  beneficial  in  areas  of  low 
threat  frequency. 

B.  RECOMMENDATIONS 

The  following  recommendations  are  offered  to  future 
researchers : 

1.  Remove  the  MAXPROB  I  and  natural  regression  strategies 
of  the  Preisendorf er  methodology  from  further  consideration 

in  the  forecasting  of  marine  horizontal  visibility. 

2.  Delete  one  of  the  two  threshold  models  evaluated 
in  this  study,  and  investigate  additional  thresholding 


techniques  based  on  the  Beta  distribution  and  the  maximum- 
likelihood-of-detection  criteria . 

3.  Revise  the  current  three-category  visibility  scheme 
to  a  two-category  scheme  where  visibility  categories  I  and 
II  are  combined.  This  should  be  particularly  beneficial  in 
those  homogeneous  areas  with  extremely  low  frequencies  of 
visibility  category  I  and  II  observations. 

4.  Expand  the  initial  set  of  potential  predictors  to 
include  air-sea  temperature  differences  as  well  as  additional 
derived  predictors  such  as  the  advections  and  gradients  of 
temperature,  vorticity  and  moisture  in  order  to  more  fully 
simulate  the  physical  processes  associated  with  poor  marine 
visibility.  Additionally,  include  TAU-00  and  TAU-24  model 
output  parameter  fields  as  potential  predictors  in  future 
evaluations  of  TAU-24  and  TAU-48  MOS  forecasts. 

5.  Evaluate  0000GMT  data  sets  to  determine  the  effect 
of  nighttime  conditions  on  both  visibility  observations  and 
NOGAPS  model  output  parameters. 

6.  Investigate  new  procedures  to  determine  the  number 
of  equally  populous  predictor  intervals.  The  following 
procedure  [Preisendorf er ,  1984]  is  proposed: 

a.  To  establish  the  number  of  equally  populous  predic¬ 
tor  intervals  for  any  predictor,  consider  a  bivariate 
predictand/predictor  [Preisendor fer ,  1983a].  Start  with 
m  =  1  and  find  the  potential  predictability  (PP)  for  the 
resultant  plot,  call  it  "PP(1)."  In  general,  PP(m)  is  the 


PP  for  the  general  case  of  m.  Successively,  find  PP(m) 
for  m  =  1,2,...,  and  continue  to  subidivde  the  predictor 
range  as  long  as  PP  increases:  PP(m)  <  PP(m+l).  Stop  at 
PP(m)  if  PP(m+l)  <_  PP(m)  or  if  PP(m+l)  <  PP(m+l|96),  where 
the  later  is  defined  by  Preisendor f er  (1983a)  and  denoted 
by  "PP(96)."  This  last  condition  avoids  sparse  bivariate  data 
plots,  caused  by  too  large  an  m.  It  was  Karl's  1984  experi¬ 
ence  that  five  to  eight  equally  populous  predictor  intervals 
are  sufficient  for  all  predictors.  Hence  m,  for  each  pre¬ 
dictor,  is  expected  to  be  in  this  neighborhood. 

b.  Order  the  set  of  available  predictors  in  descending 
value  of  potential  predictability  (PP) .  Break  ties  with 

AO  (PP  and  AO  are  defined  by  Preisendorf er  (1983a)).  AO 
is  the  actual  skill,  after  the  prediction  has  been  made. 

c.  The  first  predictor  is  that  with  the  greatest 

PP.  Compute  associated  AO  and  Al.  Call  them  AO^  and  Al^. 

7.  Associated  with  recommendation  6.  above,  improve 
the  predictor  selection  procedure  as  follows: 

a.  Suppose  k-1  predictors  have  been  chosen,  let 
them  be  X^,...,X^_^.  Let  Y  be  a  new  predictor  candidate. 

Admit  Y  as  the  kth  predictor  if  the  three  following  conditions 
are  satisfied: 

(1)  Functional  dependence  (y|X^)  <  functional 
dependence^  (Y|x^;05)  for  i  =  l,...,k-l 
(i.e.,  the  functional  dependence  of  X^ 
and  Y  is  not  significantly  large  for  each 
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i  =  l,...,k-l.  Find  functional  dependence 
(Y|xi)  and  functional  dependence  (YjX^;05) 
as  described  by  Preisendorf er  (1983c)). 


(2) 

> 

o 

i — \ 

i 

O 

< 

A 

and  Al(k) 

>  Al(k-1} 

(3) 

A0(k) 

>  AO (96) 

and  Al ^ k  ^ 

>  Al (05) 

All  three  conditions  must  hold  for  admittance  of  Y  to  the 
predictor  set. 

b.  A  less  stringent  predictor  selection  process 
would  be  to  form  functional  dependence  (Y^|x^),  where  X^, 
i  =  l,...,k-l  are  the  selected  predictors,  and  the  Y., 

j  =  1 , . . . , q  are  the  as-yet  unselected  predictors.  Here 
(k-l)+q  =  p,  the  original  number  of  potential  predictors. 
Next,  form  |min|max  functional  dependence  (Y.|x.)|j.  This 

j  i _ 3  1 

fixes  that  Y_.  for  which  functional  dependence  (Y^jx^)  is 

the  least  possible  of  the  maximum  functional  dependence 

values  over  the  present  X  set.  This  makes  the  best  out  of 

the  worst  case  of  functional  dependence  (i.e.,  select  the 

Y.  farthest  from  the  set  of  X.). 

3  i 

c.  Continue  to  repeat  step  7a.  above  until  all 
potential  predictors  are  used  up  (the  critical  values  of 
A0(96)  and  Al(05)  are  as  defined  by  Preisendorf er  (1983b). 
Another  reason  for  stopping  may  be  that  allotted  CPU  time 
used  up  before  the  predictors. 


is 


8.  Investigate  a  further  and  more  complete  verification 
of  the  forecast  arrays  and  equations  presented  in  this  study 


utilizing  available  1984  data  sets.  Specif ically ,  utilize 
the  1984  data  set  as  an  entire  independent  test  case  without 
first  removing  a  portion  of  the  data  for  use  as  a  dependent 
forecast  array/equation  training  set.  Additionally,  gener¬ 
ate  an  additional  set  of  forecast  arrays  and  equations 
based  on  combined  1983  and  1984  FATJUNE  data  and  evaluate 
the  statistical  stability  of  the  equations  as  different 
years  of  data  are  merged  into  a  larger  data  base. 


APPENDIX  A 


LINEAR  REGRESSION  AND  THRESHOLD  MODELS 

A.  LINEAR  REGRESSION 

The  linear  regression  techniques  used  in  this  study 
expand  upon  and  slightly  modify  those  first  presented  by 
Karl  (1984).  In  this  study,  two  separate  least-squares, 
multiple  linear  regression  software  programs;  referred  to 
as  the  BMDP2R — Stepwise  Regression  and  BMDP9R — All  Possible 
Subsets  Regression  computer  programs  in  the  BMDP  Statistical 
Software  [University  of  California,  1983],  were  used. 

The  independent  variable  selection  procedure  employed  in 
the  BMDP2R  program  is  referred  to  as  a  forward,  step-wise 
selection  process  where  predictors  are  selected  from  a  large 
group  of  available  potential  independent  variables  based 
upon  the  highest  correlation  with  the  dependent  predictand 
(visibility  in  this  study) .  This  correlation  is  calculated 
based  upon  certain  F-to-enter  and  F-to-remove  limits,  where 
F  is  a  ratio  which  tests  the  significance  of  the  coefficients 
of  the  predictors  in  the  regression  equation. 

The  regression  model  fitted  to  the  data  is 

y  =  a  +  b.x.  +  b0x_  +  ...  +  b  x  +  e 
112  2  n  n 


where ; 


y  =  the  dependent  variable  (predictand)  which 
can  be  either  a  continuous  function  or  a 
discrete  value 

x^,...,xn  =  the  independent  variables  (predictors) 

b^,...,b  =  the  regression  coefficients 

a  =  the  intercept 

p  =  the  number  of  independent  variables 
e  =  the  error  with  mean  zero 

The  predicted  value  y,  and  the  general  form  of  the  resulting 
equation,  is 

A 

y  =  a  +  blXl  +  b2x2  +  . . .  +  bnxn 

The  step-wise  selection  of  predictors  continues  until  there 

are  no  predictors  remaining  which  meet  the  requisite  F-to- 

enter  criteria.  The  regression  equation  generated  by  the 

BMDP2R  program  is  outputted  at  each  regression  step  where 

variables  are  selected  as  independent  predictors,  along  with 

its  corresponding  R  value  (the  correlation  of  dependent 

~  2 

variable  y  with  the  predicted  value  y)  and  R  value.  The 

resulting  equation  sets  are  reviewed,  and  that  equation  con- 

2 

taining  only  those  predictors  which  increased  R  by  at  least 
0 . 01  are  retained  for  application. 

The  procedure  employed  with  the  BMDP9R  program  varies 
from  that  of  the  BMDP2R,  in  that  a  "best"  possible  subset, 
derived  independently  of  variable  or  variable  sequence,  is 
calculated  from  the  group  of  potential  predictors.  Once  this 


"best"  subset  is  identified,  a  linear  regression  equation  is 
fitted  to  the  data,  based  only  upon  those  selected  predic¬ 
tors,  in  a  fashion  identical  to  that  for  the  BMDP2R  program. 
The  "best"  possible  subset  is  calculated  by  a  Furnville- 
Wilson  algorithm  which  provides  the  user  with  a  variety  of 
subordinate  subsets  in  addition  to  the  "best"  subset.  Three 
criteria  are  available  to  define  the  "best"  possible  subset 

as  a  function  of  independent  variables  (predictors)  and  a 

2 

dependent  variable  (predictand) :  the  sample  R  ,  the  adjusted 
2 

R  and  Mallow's  Cp.  For  this  study,  the  Mallow's  Cp  criteria 
is  defined  as: 

Cp  =  RSS/S  -  (N  -  2P '  ) 

where : 

RSS  =  the  residual  sum  of  the  squares  for  the 
new  subset  being  tested 

S  =  the  residual  mean  square  based  on  the 

linear  regression  using  all  independent 
variables 

P'  =  the  number  of  variables  in  each  subset 
N  =  the  total  number  of  cases 

For  this  method,  "best"  is  defined  as  the  smallest  Cp  value. 

Independent  variable  selection  for  the  BMDP9R  program 
begins  with  a  general  screening  of  the  entire  set  of  potentia 
predictors.  Variables  which  are  identified  as  redundant. 
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linear  combinations  of  other  variables,  with  respect  to  the 
predictand,  in  this  general  screening  are  deleted  from  further 
consideration.  The  t  statistics  for  the  coefficients  which 
minimize  the  Cp  value  for  each  reviewed  subset  identifies 
the  "best"  subset.  The  number  of  predictors  assigned  to  each 
subset  can  be  predefined  and  for  this  study  each  subset 
equation  was  required  to  have  six  predictors. 

The  role  of  regression,  once  appropriate  predictor  varia¬ 
bles  have  been  selected,  is  simply  that  of  dimension  reduction 
(representing  a  multivariate  structure  by  a  univariate 
proxy  which  constitutes  a  classif icatory  or  predictive 
index) .  This  proxy  takes  the  form  of  a  polynomial,  linear 
in  its  coefficients,  of  the  components  of  the  multivariate 
structure.  The  problem  now  becomes  one  of  determining  the 
form  of  the  state  conditional  distributions  (one  for  each 
group  of  interest;  e.g.,  one,  two  and  three  for  visibility 
categories  I,  II  and  III,  as  used  in  this  study).  Once  an 
appropriate  form  has  been  selected,  it  remains,  then,  to 
determine  the  parameters  of  the  class  conditional  distribu¬ 
tions  (e.g.,  means  and  variances)  and  then  apply  an  appro¬ 
priate  decision  criterion  or  threshold  model. 

B.  THRESHOLDS  [Lowe,  1984a] 

1 .  Notation 

E  =  an  event;  this  is  an  indicator  variable  which 
when  E  =  1,  the  threatening  event  occurs,  and 
when  E  =  0,  the  non- threatening  event  occurs. 
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C  =  the  classification  of  an  unknown  event  wh-'ch 
when  C  =  1,  the  event  is  classified  as  a 
threat,  and  when  C  =  0,  the  event  is  classi¬ 
fied  as  a  non- threat. 

P [ E  = 1 ]  =  unconditional  probability  of  occurrence  of 

threat . 

P[E=0]  =  unconditional  probability  of  occurrence  of 

non- threat . 

Error  of  the  1st  kind  (false  alarm)  [C  =1  n  E  =  0]  . 

Error  of  the  2nd  kind  (miss)  [C  =  0  n  E  =  1]  . 

P[C=lnE=0]  =  joint  probability  of  an  error  of  the  1st 
kind . 

P[C=0nE=l]  =  joint  probability  of  an  error  of  the 
2nd  kind. 

P [C  = 1 1 E  =  0]  =  class  conditional  probability  of  misclassi- 

fying  a  non- threat. 

P[C=0|E=1]  =  class  conditional  probability  of  misclassi- 

fying  a  threat. 

P[C  =1  n  E  =  0]  -  P  [C  =  1 1  E  =  0  ]  P  [E  =  0  ]  . 

P  [C  =  0  n  E  =  1]  =  P  [C  =  0  |  E  =  1]  P  [E  =0]  . 

z  =  a  value  of  the  predictive  index  (equivalent 
to  y ,  above) . 

Z  =  range  of  the  predictive  index  on  the  real  line. 


For  a  dichotomous  problem,  Z  is  into  two  parts  ZQ ,  Z^, 


C  =  0  if  z  s  Zg 


C  =  1  if  z  e  Z1 


The  decision  regions  are  mutually  exclusive  and  exhaustive 
(i.e.,  ZQ  =  0  and  Z  =  Zq  u  Z^)  . 
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Thresholds  =  boundary (s)  between  decision  regions. 


p(z  |e  =0) 


class  conditional  density  of  z  given 
that  E  =  0. 


p(  z  |  E  =  1] 


class  conditional  density  of  z  given 
that  E  =  1. 


A(z)  =  p(z|E=l)/p(z|E=0)  =  the  maximum  likelihood 


ratio  (i.e.,  the  ratio  of  class  conditional 
densities ) . 


Pe  =  p{  [C  =  1  n  E  =0]  u  [C  =0  n  E  =  1]  }  =  the  total 


probability  of  error. 


2.  Minimum  Probability  of  Error  Criterion 


Pe  =  probability  of  an  incorrect  classification. 


Pe  =  p  [C  =  1 1  E  =  0  ]  p  [E  =  0  ]  +  p  [C  =  0  |  E  =  1]  p[E=l] 


where  p[E  =1]  +  p[E=0]  =  1.  NOte  that  the  events  E  =  1 


and  E  =  0  are  mutually  exclusive  and  exhaustive.  The  objec¬ 


tive  is  to  select  decision  regions  (thresholds)  so  as  to 


minimize  p 


p[C=0|E=l]  =  /  p(z|E=l)dz  =  the  probability  of 

zeZ 


misclassif ying  E  =  1. 


p  [C  =  0  |  E  =  1  ]  =  /  p  (  z  |  E  =  1 )  dz  +  /  p  (  z  |  E  =  1)  dz 


/  p ( z ] E  =  1 ) dz 
zcZ. 


p  (C  =  0  |  E  =  1]  = 


1-  /  p  (  z  |  E  =  1 )  dz 

Zf.  Z, 


p(C  =  1  |E  =0]  = 


/  p(z|E=0)dz 
z<=.  Z. 


these  are 
substituted 
into  the 
expression 

for  pQ 


'  .o: 
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then , 


p  =  p  [E  =  0  ]  /  p(z|E=0)dz  +  p  [E  =1)  [  1  -  /  p(z|E  =  l)dz] 


ze  Z . 


zeZ. 


and  algebraic  rearrangement  yields, 


p  =  p[E  =1]  -  /  { p  [  E  =  0  ]  p  (  z  |  E  =  0 )  -  p  [E  =  1)  p  (  z  |  E  =  1)  dz} 

zeZ. 


In  order  to  minimize  p^,  Z^  (the  decision  region  for  C  =  1) 
will  include  all  those  values  of  z  for  which  the  integrand 
in  the  expression  for  pg  will  be  negative.  The  decision 
regions  can  be  symbolically  represented  as  follows: 


Z_  =  {z:  p[E=0]  p(z|E=0)  -  p  [E  =  1]  p(z|E=l)  >  0} 


Z1  =  {z:  p  [E  =  0]  p(z|E=0)  -  p[E=lJ  p(z|E=l)  <  0} 


An  alternative  representation  is  given  by. 


Zq  =  {z:  p  [E  =  0  ]  p(z|E=0)  >  p[E=l]  p(z|E=l)} 


=  (z:  p  [E  = 0]/p[E  = 1]  >  p(z|E=l)/p(z|E=0)} 


Likewise, 


Z1  =  (z:  p[E  =  0]/p[E  =  1]  <  p  (  z  |  E  =  1 )  /p  (  z  |  E  =  0)  } 
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These  statements  can  be  combined  to  give, 


c=l 

p  (  z  |  E  =  1) /p  (  z  |  E  =  0  )  =  A  (  z )  *  p[E  =0]/p[E  =  1] 

c-0 

Thresholds  are  the  value (s)  of  z  for  which 


A  ( z)  =  p[E  =0]/p[E  =  1] 


This  equation  can  be  solved  for  z  either  analytically  or 
numerically  depending  on  the  forms  of  the  density  functions 
3 .  Threshold  Cases 

In  order  to  examplify  the  model,  the  assumption  is 

made  that  the  class  conditional  distributions  are  Gaussian. 

There  are  essentially  three  distinct  cases  that  can  arise. 

a.  Case  I:  Equal  variances;  different  means 

(Referred  to  as  the  equal  variance  model  (EVAR) 
in  the  text) 

p  (  z  |  E  =  1)  =  k  exp{  (-1/2)  (z  -Ul)2/o2} 
p ( z  |  E  =  0 )  =  k  exp{ (-1/2)  (z  -  pQ) 2/o2} 


where : 


k  = 


, ~l/2  -1 
(  2  TT )  0 


exp{  (-1/2)  (z  -u1)2/a2}  C>1  pQ 
exp{ (-1/2) (z  -Un)2/o2}  '  P1 


Note  that  in  this  situation  there  are  two  thresholds.  The 


group  having  the  smaller  variance  will  lie  between  the  two 
thresholds . 


The  thresholds  shown  are  typical  of  a  situation  where  p^  <  p^ 
Note  that  these  thresholds  lie  between  the  two  intersections 
of  the  densities.  If  the  inequality  of  prior  probabilities 
were  reversed,  the  thresholds  would  lie  outside  of  the 
region  between  the  two  density  intersections.  Further  note 
that  the  decision  region  for  the  group  having  the  lesser 
variance  lies  between  the  thresholds. 

c.  Case  III:  General  Solution  (Referred  to  as  the 
Quadratic  Model  (QUAD)  in  the  text) 

p  ( z  |  E  =  1 )  =  k/o1  exp{  (-1/2)  (z  -  y:)  2/o21) 


p ( z | E  =  0) 


k/o  exp{ (-1/2)  ( z  -  v  )  2/o2  } 


A  (z) 


exp{ 1/2  (— — 

°0 


where  k 


(2  tt) 


-1/2 


Algebraic  manipulation  produces 


2  2 
(o^-c^) 


2  (o0ul 


2  . 
~°lV 


c 

2  2  2  2  2  2 

+  t(a1M0-o0u1)  ~  2a()01  In  (p0a1/p1aQ)  ] 

c 


which  is  recognizable  as  a  quadratic  equation  in  z. 
z*  =  -b  ±  (b2  -  4ac)1/2/2a 

where : 


b  -  2(0^  -  ojy0) 

=  (o2y2  -  o2y2)  -  2o2y2  In  (pno,/Plon) 


c 


The  remarks  given  for  the  figures  in  cases  I  and  II  are  also 
applicable  here.  More  often  than  not,  only  one  of  a  pair  of 
thresholds  induced  by  differing  variances  will  be  of  real 
interest.  If  the  variances  of  the  two  groups  are  radically 
different,  then  both  members  of  the  threshold  pair  become 
important . 

4 .  The-Maximum- Likelihood- of -Detect ion  Criteria 

For  this  specific  model  the  following  background  is 
provided : 


event  space:  2'  mutually  exclusive  populations 

tTq,  tt^  forecast  decision  space:  2  possible  forecasts 


dg  is  a  correct  forecast  if  actually  occurs 
d-^  is  a  correct  forecast  if  actually  occurs 

Problem:  select  the  decision  rule  d(z)  which  maps 

the  observation  space  Z  into  some  forecast  space 
in  some  optimal  manner. 
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Z  may  be  an  observed  variable  or  it  may  be  an 
univariate  index  derived  from  a  number  of  variables. 

For  this  two  decision  problem,  Z  is  partitioned 
into  two  parts,  Zq  and  Z^. 

d ( z )  =  dQ  if  z  c  ZQ 
d(z)  =  d1  if  z  £  Zj 

where  Z  n  Z  =  0  and  Zn  u  Zn  =  Z 
0  1  0  1 

The  maximum-likelihood-of-detection  criteria  repre¬ 
sents  the  simplest  decision  model.  The  basic  involves 
selecting  the  forecast  (decision)  corresponding  to  the  obser¬ 
vation  (signal)  which  is  the  most  likely  symptom  of  the  event 
subsequently  observed.  Consider  the  following  example: 

problem:  diagnose  disease  A  or  disease  B. 

The  observed  symptoms  occur  with  probability  0.75 
for  A  and  0.1  for  B.  By  the  maximum-likelihood-of-detection 
criteria  (MLDC) ,  diagnose  disease  A  because  A  is  the  most 
likely  cause  of  the  observed  symptoms  (if  there  is  no  more 
information) .  But  if  we  know  that  A  is  rare  and  B  is  common, 
the  above  decision  may  not  be  optimal  and  MLDC  may  not  be 
appropriate.  MLDC  requires  only  that  we  know  the  event 
conditional  probability  density  functions  of  the  observations. 
That  is  : 
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"AD-A1S1  935  HN  EVALUATION  OF  DISCRETIZED  CONDITIONAL  PROBABILITY 
AND  LINEAR  REGRESSIO.  .  <U)  NAVAL  POSTGRADUATE  SCHOOL 
HONTEREV  CA  N  DIUNIZIO  SEP  84 


UNCLASSIFIED 


F/G  12/1 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

national  miRfAo  of  c  unhappy  iq**  a 


decision  rule:  d{z) 


|  d1  if  p  (  Z  |  TT  1  )  >  p  (  Z  |  TT  Q  ) 
|  dg  if  p  (  Z  |  TT  1 )  <  p  (  Z  j  TT  0  ) 

In  the  following  development  the  Gaussian  density 
is  used  to  exemplify  the  model. 

2—2 

p(z|TTn)  =  1//21H77  exp{-l/2  (— — ^-)2} 

u  U  Oq 

Z  —  2 

p  ( z  l  tt  ^ )  =  l/t^TTO^  exp{-l/2(— - — -)  2  } 

P  (  Z  I  TT  } 

definition:  likelihood  ratio  A(z)  =  , - r 

P(Z  |  TTq) 

for  convention  sake  we  assume  z^  >  Zg, 


o 


2 

1 


> 


o 


2 

0 


*  note  the  class  having  the  largest  variance  has  a 
bifurcated  decision  region. 


In  the  case  where  the  variances  are  equal,  the 
situation  simplifies  considerably. 


2°2(zi  ‘  V  +°2(zo 


2  —  —  2—2 
2oZ(z1  -  z 0 )  -a  (zZ 


2z 


(Z1  +  V 


0 


0 


(z,  +  zn) 


It  is  obvious  that  z*  is  simply  the  average  of  the 
means  of  the  class-conditional  distributions  and  is  found 
at  the  intersections  of  the  two  density  curves. 

In  the  foregoing,  normal  class  conditional  distribu¬ 
tions  were  assumed.  This  was  done  because  the  Gaussian  form 
admits  of  a  rather  clean  analytical  solution.  However,  the 
general  concept  of  the  minimum  probable  error  decision 
criteria  may  be  applied  to  any  form  of  density  function. 
Indeed,  the  density  function  of  one  group  need  not  even  be 
the  same  form  as  that  for  another  group  (one  might  be  exponen 
tial  and  the  other  Gaussian) .  The  difficulty  with  most  no n- 
Gaussian  forms  is  that  they  seldom  admit  of  closed  analytical 
forms  and  require  numerical  means  in  determination  of 


thresholds . 


APPENDIX  B 


NOGAPS  PREDICTOR  PARAMETERS  AVAILABLE  FOR 
NORTH  ATLANTIC  OCEAN  EXPERIMENTS 

Area:  Entire  North  Atlantic  Ocean  and  Mediterranean  Sea 


Model  output 

time:  1200GMT  (TAU-00) 

15  May--7  July  1983 

Model  output 

Descriptive  name  of  parameter 

parameter 

DlOOO 

1000  mb  geopotential  height 

D925 

925  mb  geopotential  height 

D850 

850  mb  geopotential  height 

D700  * 

700  mb  geopotential  height 

D500 

500  mb  geopotential  height 

D400  * 

400  mb  geopotential  height 

D300  * 

300  mb  geopotential  height 

D250  * 

250  mb  geopotential  height 

TAIR 

Surface  air  temperature 

T1000 

1000  mb  temperature 

T9  2  5 

925  mb  temperature 

T700  * 

700  mb  temperature 

T500 

500  mb  temperature 

T400  * 

400  mb  temperature 

T300  * 

300  mb  temperature 

T250  * 

250  mb  temperature 

EAIR 

Surface  vapor  pressure 

E1000 

1000  mb  vapor  pressure 

E925 

925  mb  vapor  pressure 

E850 

850  mb  vapor  pressure 

E700  * 

700  mb  vapor  pressure 

E500 

500  mb  vapor  pressure 

UBLW 

Boundary  layer  zonal  wind  component 

99 


U1000 
U925 
U850 
U700  * 
U500 
U400  * 
U300  * 
U250  * 
VBLW 

V1000 

V925 

V850 

V700  * 

V500 

V400  * 

V300  * 

V250  * 

VOR925  ** 

VOR500  ** 

PS 

SMF 

PBLD 

STRTFQ 

STRTTH 

SHF 

ENTRN 

DRAG  ** 


1000  mb  zonal  wind  component 

925  mb  zonal  wind  component 

850  mb  zonal  wind  component 

700  mb  zonal  wind  component 

500  mb  zonal  wind  component 

400  mb  zonal  wind  component 

300  mb  zonal  wind  component 

250  mb  zonal  wind  component 

Boundary  layer  meridional  wind 
component 

1000  mb  meridional  wind  component 

925  mb  meridional  wind  component 

850  mb  meridional  wind  component 

700  mb  meridional  wind  component 

500  mb  meridional  wind  component 

400  mb  meridional  wind  component 

300  mb  meridional  wind  component 

250  mb  meridional  wind  component 

925  mb  vorticity 

500  mb  vorticity 

Surface  pressure 

Surface  moisture  flux 

Planetary  boundary-layer  depth 

Percent  stratus  frequency 

Stratus  thickness 

Surface  heat  flux 

Entrainment  at  top  of  marine 
boundary- layer 

Drag  coefficient  (CD) 


Derived  parameters 

DTDP  Vertical  gradient  of  temperature 

(1000-925  mbs) 

DEDP  Vertical  gradient  of  vapor  pressure 

(1000-850  mbs) 

DUDP  Vertical  gradient  of  zonal  wind 

(1000-850  mbs) 


100 


DVDP 


RH 

TV 

DDVDP 

DVRTDP 

ESUM 

EPRD 


Vertical  gradient  of  meridional  wind 
(1000-850  mbs) 

Surface  relative  humidity 

Virtual  temperature 

Vertical  gradient  of  geopotential 
height  (1000-850  mbs) 

Vertical  gradient  of  vorticity 
(500-925  mbs) 

Sum  of  vapor  pressures 
(1000-850  mbs) 

Product  of  vapor  pressures 
(1000-850  mbs) 


Area:  Entire  North  Atlantic  Ocean  and  Mediterranean  Sea 

Model  forecast  projection:  1200GMT  (TAU-24) 

15  May — 7  July  1983 


Model  output 
parameter 

Descriptive  name  of 

parameter 

D1000 

1000  mb  geopotential  height 

D925 

925  mb  geopotential 

height 

D850 

800  mb  geopotential 

height 

D700  * 

700  mb  geopotential 

height 

D500 

500  mb  geopotential 

height 

D400  * 

400  mb  geopotential 

height 

D300  * 

300  mb  geopotential 

height 

D250  * 

250  mb  geopotential 

height 

TAIR 

Surface  air  temperature 

T1000 

1000  mb  temperature 

T92  5 

925  mb  temperature 

T700  * 

700  mb  temperature 

T500 

500  mb  temperature 

T400  * 

400  mb  temperature 

T300  * 

300  mb  temperature 

T250  * 

250  mb  temperature 

EAIR 

Surface  vapor  pressure 

ElOOO 

E925 

E850 

E700  * 

E500 

UBLW 

U1000 

U925 

U700  * 

U500 

U400  * 

U300  * 

U250  * 

VBLW 

VIOOO 

V925 

V850 

V700  * 

V500 

V400  * 

V300  * 

V250  * 

VOR925 

VOR500 

PS 

SMF 

PBLD 

STRTFQ 

STRTTH 

SHF 

ENTRN 


1000  mb  vapor  pressure 

925  mb  vapor  pressure 

850  mb  vapor  pressure 

700  mb  vapor  pressure 

500  mb  vapor  pressure 

Boundary  layer  zonal  wind  component 

1000  mb  zonal  wind  component 

925  mb  zonal  wind  component 

700  mb  zonal  wind  component 

500  mb  zonal  wind  component 

400  mb  zonal  wind  component 

300  mb  zonal  wind  component 

250  mb  zonal  wind  component 

Boundary  layer  meridional  wind 
component 

1000  mb  meridional  wind  component 

925  mb  meridional  wind  component 

850  mb  meridional  wind  component 

700  mb  meridional  wind  component 

500  mb  meridional  wind  component 

400  mb  meridional  wind  component 

300  mb  meridional  wind  component 

250  mb  meridional  wind  component 

925  mb  vorticity 

500  mb  vorticity 

Surface  pressure 

Surface  moisture  flux 

Planetary  boundary-layer  depth 

Percent  stratus  frequency 

Stratus  thickness 

Surface  heat  flux 

Entrainment  at  top  of  marine 
boundary- layer 

Drag  coefficient  (CQ) 

Total  amount  (mm.)  of  model 
precipitation  in  the  last  six  hours 


DRAG 

PRECIP 


SHWRS 

INSTAB 

DIV925 


Total  amount  (mm.)  of  model  precipita¬ 
tion  associated  with  cumulus  convection 
in  the  last  six  hours 

Boundary  layer  inversion  instability 

925  mb  Divergence 


B.  Derived  parameters 


DTDP 

DEDP 

DUDP 

DVDP 

RH 

TV 

DDVDP 

DVRTDP 

ESUM 

EPRD 


Vertical  gradient  of  temperature 
(1000-925  mbs) 

Vertical  gradient  of  vapor  pressure 
(1000-850  mbs) 

Vertical  gradient  of  zonal  wind 
(1000-850  mbs) 

Vertical  gradient  of  meridional  wind 
(1000-850  mbs) 

Surface  relative  humidity 

Virtual  temperature 

Vertical  gradient  of  geopotential  height 
(1000-850  mbs) 

Vertical  gradient  of  vorticity 
(500-925  mbs) 

Sum  of  vapor  pressures 
(1000-850  mbs) 

Product  of  vapor  pressures 
(1000-850  mbs) 


III.  Area:  Entire  North  Atlantic  Ocean  and  Mediterranean  Sea 


Model  forecast  projection: 


1200GMT  (TAU-48) 

15  May--9  July  1983 


A. 


Model  output 
parameter 

D10  00 
D92  5 
D850 
D700  * 

D500 
D400  * 


Descriptive  name  of  parameter 

1000  mb  geopotential  height 
925  mb  geopotential  height 
850  mb  geopotential  height 
700  mb  geopotential  height 
500  mb  geopotential  height 
400  mb  geopotential  height 
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'0 


j 


D300 

* 

300  mb  geopotential  height 

D250 

* 

250  mb  geopotential  height 

TAIR 

Surface  air  temperature 

T1000 

1000  mb  temperature 

T9  2  5 

925  mb  temperature 

T700 

* 

700  mb  temperature 

T500 

500  mb  temperature 

T400 

* 

400  mb  temperature 

T300 

* 

300  mb  temperature 

T250 

* 

250  mb  temperature 

EAIR 

Surface  vapor  pressure 

E1000 

1000  mb  vapor  pressure 

E9  25 

925  mb  vapor  pressure 

E850 

850  mb  vapor  pressure 

E700 

★ 

700  mb  vapor  pressure 

E500 

500  mb  vapor  pressure 

UBLW 

Boundary  layer  zonal  wind  component 

UlOOO 

1000  mb  zonal  wind  component 

U925 

925  mb  zonal  wind  component 

U850 

850  mb  zonal  wind  component 

U700 

* 

700  mb  zonal  wind  component 

U500 

500  mb  zonal  wind  component 

U400 

* 

400  mb  zonal  wind  component 

U300 

* 

300  mb  zonal  wind  component 

U250 

* 

250  mb  zonal  wind  component 

VBLW 

Boundary  layer  meridional  wind 
component 

V1000 

1000  mb  meridional  wind  component 

V925 

925  mb  meridional  wind  component 

V850 

850  mb  meridional  wind  component 

V700 

★ 

700  mb  meridional  wind  component 

V500 

500  mb  meridional  wind  component 

V400 

* 

400  mb  meridional  wind  component 

V300 

* 

300  mb  meridional  wind  component 

V250 

★ 

250  mb  meridional  wind  component 

VOR  925 

VOR500 

PS 

SMF 

PBLD 

STRTFQ 

STRTTH 

SHF 

ENTRN 

DRAG 

PRECIP 

SHWRS 

INSTAB 

DIV925 

Derived 

DTDP 

DEDP 

DUDP 

DVDP 

RH 

TV 

DDVDP 

DVRTDP 

ESUM 


925  mb  vorticity 

500  mb  vorticity 

Surface  pressure 

Surface  moisture  flux 

Planetary  boundary-layer  depth 

Percent  stratus  frequency 

Stratus  thickness 

Surface  heat  flux 

Entrainment  at  top  of  marine 
boundary- layer 

Drag  coefficient  (CD) 

Total  amount  (mm. )  of  model  precipitation 
in  the  last  six  hours 

Total  amount  (mm.)  of  model  precipitation 
associated  with  cumulus  convection  in 
the  last  six  hours 

Boundary  layer  inversion  instability 
925  mb  Divergence 

parameters 

Vertical  gradient  of  temperature 
(1000-925  mbs) 

Vertical  gradient  of  vapor  pressure 
(1000-850  mbs) 

Vertical  gradient  of  zonal  wind 
(1000-850  mbs) 

Vertical  gradient  of  meridional  wind 
(1000-850  mbs) 

Surface  relative  humidity 

Virtual  temperature 

Vertical  gradient  of  geopotential  height 
(1000-850  mbs) 

Vertical  gradient  of  vorticity 
(500-925  mbs) 

Sum  of  vapor  pressures 
(1000-850  mbs) 

Product  of  vapor  pressures 
(1000-850  mbs) 


EPRD 


used  due  to  their 
cally  unimportant 
ibility . 

used  due  to  loss  o 
transfer  from  tape 
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APPENDIX  C 


SKILL  AND  THREAT  SCORES,  DEFINITIONS  [Karl,  1984 


S  T 


V  W 


V  z 


1  2  3 

OBSERVED 


Total  =  R+S  +  T  +  U  +  V  +  W  +  X  +  Y+  Z 
PI  =  (R+U+X) /Total  P 3  =  (T+W+Z ) /Total 


P2  =  (S+V+Y) /Total 


PN  =  greatest  of  Pi,  P2  or  P3 


Raw  Scores 


AO  =  %  correct  =  ( X+V+T) /Total 


Al  =  one-class  error  =  (U+S+Y+W) /Total 

TSl  =  Threat  score  for  visibility  category  I 
=  X/ (R+U+X+Y+Z) 

TS2  =  Threat  score  for  visibility  category  II 
=  V/ (S+V+Y+U+W) 

T312  =  Threat  score  for  visibility  categories  I  and  II 


(X+V) / (Total-T) 
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TS12  is  designed  to  represent  the  skill  of  forecasting 
visibility  categories  I  and  II  as  separate  categories, 
rather  than  their  skill  as  a  combined  category,  which 
would  be  (U+V+X+Y) / (Total-T) . 


Adjusted  scores 

AAO  =  ( AO-PN) / ( 1-PN) 

ATS1  =  (TS1-P1) /(1-PI) 

ATS 2  =  (TS2-P2) /(1-P2) 

=  (TS12- ( P1+P2 ) )/(l-(Pl+P2) ) 


ATS  12 


APPENDIX  D 


BMDP  LINEAR  REGRESSION  EQUATION  PREDICTOR  SETS, 

NORTH  ATLANTIC  OCEAN  EXPERIMENTS  (PR+BMD  MODEL) 

I.  Area  2,  TAU-00  (BMDP  P2R) 

BMDl  =  2.842  -  0.21767*E850  +  0 . 837882E-05*D500 

+  0 . 0 329 3*SHF  +  7.057*DTDP  +  0.05872*ESUM  ** 

BMD2  =  -10.4469  +  0.11854*EAIR  +  0.10124*SMF 

-  0 . 07409*T9 2 5  -  0.16481*E925 

BMD3  =  3.47713  -  0.22482*EM  +  45 . 06116*DEDP 

+  0 . 00521*EPRD 

II.  Area  2,  TAU-24  (BMDP  P9R) 

BMDl  =  -20.9733  -  0.20905*E850  -  0 . 078694*T925 

+  0 . 0533674*SHF  -  0 . 0316725*INSTAB 
+  0 . 0862939*TV  +  0 . 08629 83*ESUM 

BMD2  =  2.68106  +  0.0356103*TM  +  0 . 53048E-04*V500 

-  0. 141302*E925  +  22.0764*DEDP 

+  0 . 012561 8 *DDVDP  +  0 . 00563327*EPRD  ** 

BMD3  =  -35.2882  +  0.0381891*PS  +  0 .0273575*T500 

+  0 .00449538*PBLD  -  0 . 00625203*STRTFQ 

-  0 . 0083686*STRTTH  +  0 . 00272894*DTDP 

III.  Area  2,  TAU-48  (BMDP  P9R) 

BMDl  =  -37.7157  -  0 . 14 70 84 *E500  -  0 . 0 89 7567*T9 25 

-  0 . 128407*E925  +  0.022881*SHF  +  0 . 00860574*RH 
+  0 . 145537*TV 
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BMD2  =  1.85487  +  0.0777253*TM  -  0 . 0266753*E850 

-  0 .0000390116*U500  -  0 . 0000366663*V500 

+  0 .0240246*DDVDP  +  0 . 105648*ESUM  ** 

BMD3  =  -13.9637  +  0.0160572*PS  +  0 . 00308705*PBLD 

-  0 . 0031323*STRTFQ  -  0 . 00846443*DTDP 
+  25 . 6871*DEDP  -  0 . 00296342*EPRD 

Area  3W,  TAU-24  (BMDP  P2R) 

BMDl  =  2.673  -  0.09363*E850  -  0.05101*T925 

-  0 . 20451*E925  +  0.0305*SHF  +  0.15111*ESUM  ** 

BMD2  =  1.15536  +  0.16326*EAIR  +  0.01509*SMF 

+  0 . 13014E-04*DM  +  8.08795*DEDP 
+  0 .02091*DDVDP  -  0.00788*EPRD 

BMD3  =  -18.55031  +  0.02089*PS  -  0 . 30643E-03*VBLW 

-  0 . 01151*STRTFQ  +  0 . 02772*STRTTH 

-  0 . 05736*DUDP 

Area  3W,  TAU-48  (BMDP  P9R) 

BDM1  =  1.92874  -  0 . 0719817*T925  -  0 . 2016 63 *E9 25 

+  0 . 0376905*SHF  +  7.66796*DEDP  +  0 . 182705*ESUM 

-  0. 00585094* EPRD 

BMD2  =  -33.2574  -  0.1459*E850  -  0 . 000205441*V925 

+  0 . 325802*SHWRS  +  0.0168064*RH  +  0.124434*TV 
+  0 . 0247277*DDVDP  ** 

BMD3  =  -10.1316  +  0 . 0126085*PS  -  0 . 00032403*VM 

+  0. 000112099 *U500  -  0 . 00880168*STRTFQ 
+  0.0159356*STRTTH  -  0 . 00174911*DRAG 


VI.  Area  4,  TAU-00  (BMDP  P2R) 


BMD1  = 

2.86828  +  0.4877E-04*U500 

-  0.64632E-04*V500  -  0 . 30475E-02*STRTFQ 

BMD2  = 

2.70156  +  0 . 14946E-04*DM  -  0 . 00904 *STRTTH 

+  0 . 01888*SHF  +  4 . 09  377*DTDP  ** 

BMD3  = 

2.84881  +  0. 24549 E- 03  *VBLW  -  0.10113*E850 

-  0.25666E-3*V850  +  0.03273*ESUM 

VII.  Area  4 

,  TAU-24  (BMDP  P9R) 

BMD1  = 

3.00017  +  0 . 0773367*TM  -  0 . 0000464491*V500 

-  0.103205*T925  -  3 . 7 62 67 *VRT925 

-  0.000477853*DTDP  +  9.51302*DEDP 

BMD2  = 

3.05949  -  0 . 088302*E500  +  0 . 00372204*PBLD 

-  0 .00492842*STRTFQ  +  0 . 0154289 *SHF 

-  1 . 89105  *VRT500  +  0 . 0143745*DVDP  ** 

BMD3  = 

-24.6366  +  0 . 00275966*PS  -  0 . 0549077*E850 

+  0 . 195977*T500  +  0 . 0140852*INSTAB 

+  0 . 0886378*TV  +  0 . 0231264*DDVDP 

VIII.  Area 

4,  TAU-48  (BMDP  P9R) 

BDMl  = 

-27.0959  -  0.0732462*E850  +  0.01616*T500 

-  0 . 117787*T925  +  0 . 287098*SHWRS 

-  4 . 13253*VRT925  +  0.111855*TV  ** 

BMD2  = 

-3.1619  +  0 .00678538*PS  -  0 . 0850887*E500 

-  0 . 0000502297*U500  -  0 . 0000396501*V500 

-  0 . 00 389051*STRTFQ  -  0 . 00917554*RH 
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BMD3 


=  2.08319  +  0.0771067*TM  -  0 . 0282767*E925 

-  1. 5829*VRT925  -  0 . 00114379*DTDP 
+  13 . 20  73  *DEDP  +  0 . 0289832*DDVDP 


*★ 


Equation  selected  as  predictors  in  the  PR+BMD  model 


APPENDIX  E 


BMDP  LINEAR  REGRESSION  EQUATION  PREDICTOR  SETS 
FOR  TWO- STAGE  THRESHOLD  MODELS 


I.  Area  2,  TAU-00  Threshold  Equations  (BMDP  P2R) 


a . 


VISCAT  I  vs.  VISCAT  II+III  separation 

VIS  =  0.87475  -  0 . 11042*E850  +  0.01173*SHF 

+  2 . 61984*DTDP  +  0.03863*ESUM 


b . 


VISCAT  II  vs.  VISCAT  III  separation  (EVAR  model) 

VIS  =  1.03838  +  0 . 03668*EAIR  -  0.10423*E850 

+  0 . 01560  *SHF 


c . 


VISCAT  II  vs.  VISCAT  III  separation  (QUAD  model) 

VIS  =  1.03730  +  0 . 0 3727 *EAIR  -  0.10505*E850 

+  0 . 0 1552*SHF 


II.  Area  2,  TAU-24  Threshold  Equations  (BMDP  P9R) 

a.  VISCAT  I  vs.  VISCAT  II+III  separation 

VIS  =  0.792157  -  0 . 00404 602 *SMF  +  0.0271941*TM 

-  0 .0767329*E850  -  0 . 0385021*T925 

-  0 . 0656522 *E9 25  +  0 . 0683159*ESUM 


b.  VISCAT  II  vs.  VISCAT  III  separation  (EVAR  model) 

VIS  =  -13.9955  -  0 . 0612856*E850 

+  0 . 0478582*SHF  -  0 . 0303761 *INSTAB 
-  1 . 585*VRT9  2  5  +  0.0539318*TV 
+  0 . 0102962*DDVDP 
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c.  VISCAT  II  vs.  VISCAT  III  separation  (QUAD  model) 

VIS  =  -13.1384  -  0 . 107716*E850 

+  0. 00002923 *V 500 
+  0 . 0122899*SHF  +  0.0498382*TV 
+  0 . 0105859*DDVDP  +  0 . 02329 83*ESUM 

III.  Area  2,  TAU-48  Threshold  Equations  (BMDP  P9R) 

a.  VISCAT  I  vs.  VISCAT  II+III  separation 

VIS  =  -18.0327  -  0 . 0569 882 *E8 50 

-  0 .104458*T925  -  0 . 00116092*PBLD 
+  0.0957819  *PRECIP  -  0.0705180*TV 

-  0 .0118901*DVDP 

b.  VISCAT  II  vs.  VISCAT  III  separation  (EVAR  model) 

VIS  =  -14.4404  -  0 .0670910*E850  -  0.11931*E500 

+  0 . 0157788*SHF  -  0 . 023835*DUDP 
+  0 . 0551551*TV  +  0 . 00761201*DDVDP 

c.  VISCAT  II  vs.  VISCAT  III  separation  (QUAD  model) 

VIS  =  -56.9191  -  0 . 168467*TM  -  0 . 116816*E500 

+  0 . 0201013*SHF  +  0 . 000599  864  *DTDP 
+  0 . 210713*TV  -  0 . 04 87371 *ESUM 

IV.  Area  3W,  TAU-24  Threshold  Equations  (BMDP  P2R) 

a.  VISCAT  I  vs.  VISCAT  II+III  separation 

VIS  =  0.88319  -  0 . 04039*E850  -  0.03385*T925 

-  0 . 11313*E9 25  -  0 . 02692*DVDP 
+  0 . 084  3*ESUM 
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b.  VISCAT  II  vs.  VISCAT  III  separation  (EVAR  model) 

VIS  =  0.59606  +  0 .000012576*DM 

+  0 .0151*INSTAB  +  10. 23334*DEDP 

-  0.0019117 *EPRD 

c.  VISCAT  II  vs.  VISCAT  III  separation  (QUAD  model) 

VIS  =  0.73427  -  0 . 12863*E850  +  0 . 0015056*PBLD 

+  0 . 02073*SHF  +  0.04270*ESUM 

Area  3W,  TAU-48  Threshold  Equations  (BMDP  P9R) 

a.  VISCAT  I  vs.  VISCAT  II+III  separation 

VIS  =  0.424612  -  0 . 000099 7789 *VM 

+  0.0000403063*U500  -  0 . 0360747*T925 

-  0 . 132972*E925  +  0 . 114795*ESUM 

-  0 . 00286176*EPRD 

b.  VISCAT  II  vs.  VISCAT  III  separation  (EVAR  model) 

VIS  =  -7.19287  -  0 . 000 788 62 7 *PS 

-  0 . 0685555*E850  -  0 . 05799 88*E500 
+  0 . 0222827*SHF  +  0 . 195715*SHWRS 
+  0 . 0208835*ESUM 

c.  VISCAT  II  vs.  VISCAT  III  separation  (QUAD  model) 

VIS  =  -7.28759  -  0 . 000797890*PS 

-  0 . 701888*E850  -  0 . 0547155*E500 
+  0 . 0  226117  *SHF  +  0 . 202564*SHWRS 
+  0 . 02122 83 *ESUM 

Area  4,  TAU-00  Threshold  Equations  (BMDP  P2R) 

a.  VISCAT  I  vs.  VISCAT  II+III  separation 

VIS  =  0.98483  -  0 . 75281E-03*STRTFQ 

+  0 . 70578*DTDP 


VII.  Area  4,  TAU-24  Threshold  Equations  (BMDP  P9R) 


VISCAT  I  vs.  VISCAT  II+III  separation 
VIS  =  -3.02386  -  0 . 200326*T9  25 


-  0 ,000654475*STRTFQ  -  0 . 000295813*DRAG 

-  0.466227 * VRT9 2 5  -  0 . 00192255*RH 
+  0 . 0153547*TV 


VISCAT  II  vs.  VISCAT  III  separation  (EVAR  model) 

VIS  =  1.19369  -  0 . 0709575*E500 

-  0. 000617273 *V500  +  0 . 193681*SHWRS 

-  0.000551506* DRAG  -  2 . 9 6306*VRT9 25 

-  0 . 018004  8 *DUDP 


VISCAT  II  vs.  VISCAT  III  separation  (QUAD  model) 

VIS  =  1.36198  -  0 . 07379 82*E500 

-  0. 000061033 *V500  +  0 . 185475*SHWRS 

-  3 . 05142*VRT925  -  0 . 0184565*DUDP 

-  0 . 00185381*RH 


VIII.  Area  4,  TAU-48  Threshold  Equations  (BMDP  P9R) 


VISCAT  I  vs.  VISCAT  II+III  separation 

VIS  =  -8.24765  -  0 . 00 8 599 6 *T500 

-  0 . 0333 863*T9 25  -  0 . 0177121*E925 

-  0 . 00059629*STRTFQ  +  0 . 000128675*DTDP 
+  0.0343295  *TV 


VISCAT  II  vs.  VISCAT  III  separation  (EVAR  model] 

VIS  =  1.20788  -  0 . 026601*E850 

-  0 . 0000310346*U500  -  0 . 0000480308*V500 
+  0 . 203575*SHWRS  -  3 . 2719*VRT925 
+  4 . 10651*DEDP 


*.  ‘.VS,  *.  _ 
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TABLE  III.  A  summary  of  1200GMT  observations,  15  May--09 
July  1983,  North  Atlantic  Ocean  homogeneous 
areas  shown  in  Fig.  1:  TAU-48  forecast  projection 
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TABLE  IV. 


Number  of  observations  of  three  visibility 
categories  and  95%  confidence  intervals  for 
the  dependent  and  independent  FATJUNE  1983 
data  for  the  North  Atlantic  Ocean  homogeneous 
areas  2  and  4,  for  TAU-00,  TAU-24  and  TAU-48 
and  area  3W  for  TAU-24  and  TAU-48 


TAU  00 


Area 

_2 

Total 

VISCAT  I 

VISCAT  II 

VISCAT  III 

DEP 

1912 

190  (.099) 

214  (.112) 

1508  (.789) 

95% 

IND 
C.  I . 

955 

87  (.091) 
(.086-. 107) 

103  (.108) 

( . 099- .122) 

765  (.801) 

(  .778-. 808) 

Area 

_4 

Total 

VISCAT  I 

VISCAT  II 

VISCAT  III 

DEP 

3181 

85  (.027) 

400  (.126) 

2696  (.848) 

IND 

1590 

44  (.028) 

197  (.124) 

1349  (.848 

95% 

C.I. 

(.022-. 032) 

( .116-. 135) 

(.838-. 858) 

TAU 

24 

Area 

_2 

Total 

VISCAT  I 

VISCAT  II 

VISCAT  III 

DEP 

1760 

180  (.102) 

206  (.117) 

1374  (.781) 

IND 

879 

71  (.084) 

98  (.111) 

710  (.808) 

95% 

C.I. 

(.081-. 106) 

(.103-. 127) 

(.774-. 805) 

Area 

3W 

Total 

VISCAT  I 

VISCAT  II 

VISCAT  III 

DEP 

1415 

270  (.191) 

173  (.122) 

972  (.687) 

IND 

707 

137  (.194) 

89  (.126) 

481  (.680) 

95% 

C.I. 

(  .173-.  205) 

(.109-. 137) 

(.665-. 705) 

Area 

4 

Total 

VISCAT  I 

VISCAT  II 

VISCAT  III 

DEP 

2938 

81  (.028) 

368  (.125) 

2489  (.847) 

IND 

1469 

38  (.026) 

175  (.119) 

1256  (.855) 

95% 

C.I. 

(  .022-.  032) 

( .114-. 133) 

(  .839-  .860) 

121 


TAU  48 


TABLE  IV  (CONTINUED) 


Area 

_2 

Total 

VISCAT  I 

VISCAT  II 

VISCAT  III 

DEP 

1852 

182  (.098) 

230  (.124) 

1440  (.778) 

IND 

925 

91  (.098) 

107  (.116) 

727  (.786) 

C.I. 

(  .087-. 109) 

( .109-. 133) 

(.765-. 796) 

Area 

3W 

Total 

VISCAT  I 

VISCAT  II 

VISCAT  III 

DEP 

1487 

290  (.195) 

186  (.125) 

1011  (.680) 

IND 

743 

132  (.178) 

111  (.149) 

500  (.673) 

C.I. 

(  .173-.  205) 

( .119-. 147) 

(.658-. 697) 

Area 

_4 

Total 

VISCAT  I 

VISCAT  II 

VISCAT  III 

DEP 

3056 

109  (.036) 

406  (.133) 

2541  (.831) 

IND 

1527 

45  (.029) 

196  (.128) 

1286  (.842) 

C.I . 

(  .028-. 039) 

( .122-. 141) 

( .824-. 846) 

QUA 


TABLE  VI.  A  summary  of  1200GMT  observations,  15  May--23 
June  1984,  North  Atlantic  Ocean  homogeneous 
areas  shown  in  Fig.  1:  TAU-24  forecast  projection 


TABLE  VII.  Summary  of  the  contingency  table  statistics  results 

for  all  the  models  used  in  the  North  Atlantic  homogeneous 
areas  2  and  3W  for  the  evaluated  TAU-24  model  forecast 
projection  for  15  May--23  June  1984 
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APPENDIX  G 


FIGURES 


OPERATIONAL  U.5.  NAVY 
MOOEL  OUTPUT  STATISTICS  (AOS) 
DEVELOPMENT  SCHEDULE 


VISIBILITY 

CEILING 

CLOUD  AMOUNT 
OBSTR  TO  VIS 

SIG  WAVE  HT 
WIND  WAVES 
OCEAN  SWELL 

SEC  WIND 

SFC  1LMP 

SEC  DEWPOINT 

PREC1P  OCCUR 
PRECIP  AMT 

TRW  OCCUR 

NORTH  ATLANTIC 

1987 

1988 

*  MEDITERRANEAN 

1990 

1991 

NORTH  PACIFIC 

1988 

1989 

INDIAN 

SOUTH  ATLANTIC 

1989 

1990 

1991 

1992 

SOUTH  PACIFIC 

Figure  1.  Proposed  U.S.  Navy  Model  Output  Statistics 
(MOS)  development  schedule 
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Homogeneous  areas  for  the  North  Atlantic  Ocean 
May,  June  and  July,  from  Lowe  (1984b) 


SELECTION  OE  GROUPING  SIZE 


(independent  data)  for  area  2,  TAU-00  (PR  model,  MAXPROB  II 
strategy) .  For  this  case  a  grouping  size  of  eight  was  selected 


NUMBER  OP  PREDICTORS  VS.  SCORING  TECHNIQUES  independent  data 

TflUOO  -  RRCO  2  -  flflXPROBl  -  GROUPING  SIZE  OP  8  _ _ 


Figure  4a.  Skill  diagram  and  contingency  table  results  for  area  2,  TAU-00 
(PR  model)  for  the  MAXPROB  I  strategy.  The  contingency  table 
corresponds  to  those  scores  associated  with  the  maximum  independent 
VISCAT  I  threat  score  achieved  at  the  first  predictor  level 
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VISCAT  I  threat  score  achieved  at  the  first  predictor  level 


INDEPENDENT  DATA 
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Fiqure  4c.  Skill  diagram  and  contingency  table  results  for  area  2,  TAU  00 

(PR  model)  for  the  natural  regression  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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NUMBER  OF  PREDICTORS  VS.  SCORING  TECHNIQUES 


Figure  7a.  Skill  diagram  and  contingency  table  results  for  area  2,  TAU-00 

( PR+BMD  model)  for  the  MAXPROB  I  strategy.  The  contingency  table 
corresponds  to  those  scores  associated  with  the  maximum  independent 
VISCAT  I  threat  score  achieved  at  the  first  predictor  level 


NUMBER  Of  PREDICTORS  VS.  SCORING  TECHNIQUES  independent  data 

TAUOO  -  AREA  2  -  MAXPPOBI1  -  GROUPING  SIZE  Of  7  _ _ _ 


gure  7b.  Skill  diagram  and  contingency  table  results  for  area  2,  TAU-00 

( PR+BMD  model)  for  the  MAXPROB  II  strategy.  The  contingency  table 
corresponds  to  those  scores  associated  with  the  maximum  independent 
VISCAT  I  threat  score  achieved  at  the  first  predictor  level 
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Figure  9.  Contingency  table  results  for  the  area 
TAU-00  equal  variance  threshold  model 
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Figure  10. 


Contingency  table  results  for  the  area  2, 
TAU-00  quadratic  threshold  model 
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SELECTION  OE  GROUPING  SIZE 


AO  (dependent  data)  and  the  adjusted  VISCAT  I  threat  score 
(independent  data)  for  area  2,  TAU-24  (PR  model,  MAXPROB  II 
strategy) .  For  this  case  a  grouping  size  of  ten  was  selected 
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NUMBER  OE  PREDICTORS  VS.  SCORING  TECHNIQUES  independent  data 
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igure  12c.  Skill  diagram  and  contingency  table  results  for  area  2,  TAU-24 

(PR  model)  for  the  natural  regression  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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Figure  13.  Functional  dependence,  AO/Al  statistics  and  96%/05% 
confidence  interval  values  for  area  2,  TAU-24  (PR 
model,  MAXPROB  II  strategy) 
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independent  VISCAT  I  threat  score  achieved  at  the  first 
predictor  level 
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gure  15c.  Skill  diagram  and  contingency  table  results  for  area  2,  TAU-24 

(PR+BMD  model)  for  the  natural  regression  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 


Predictor 


Figure  16.  Functional  dependence,  A0/A1  statistics  and  96%/05 
confidence  interval  values  for  area  2,  TAU-24 
( PR+BMD  model,  MAXPROB  II  strategy) 
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Figure  17.  Contingency  table  results  for  the  area  2, 
TAU-24  equal  variance  threshold  model 
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Figure  18.  Contingency  table  results  for  the  area  2 
TAU-24  quadratic  threshold  model 
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AO  (dependent  data)  and  the  adjusted  VISCAT  I  threat  score 
(independent  data)  for  area  2,  TAU-48  (PR  model,  MAXPROB  II 
strategy) .  For  this  case  a  grouping  size  of  nine  was  selected 
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Figure  21a.  Skill  diagram  and  contingency  table  results  for  area  2,  TAU-48 
(PR  model)  for  the  MAXPROB  1  strategy.  The  contingency  table 
corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  second 
predictor  level 
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corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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Figure  21c.  Skill  diagram  and  contingency  table  results  for  area  2,  TAU-48 

(PR  model)  for  the  natural  regression  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fourth 
predictor  level 
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(independent  data)  for  area  2,  TAU-48  (PR+BMD  model,  M/  PROB  II 
strategy) .  For  this  case  a  grouping  size  of  nine  was  selected 
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igure  23a.  Skill  diagram  and  contingency  table  results  for  area  2,  TAU-48 
(PR+BMD  model)  for  the  MAXPROB  I  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  first 
predictor  level 
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(PR+BMD  model)  for  the  MAXPROB  II  strategy.  The  contiri' 
table  corresponds  to  those  scores  associated  with  the  m. 
independent  VISCAT  I  threat  score  achieved  at  the  first 
predictor  level 
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Figure  23c.  Skill  diagram  and  contingency  table  results  for  area  2, 
(PR+BMD  model)  for  the  natural  regression  strategy.  The 
contingency  table  corresponds  to  those  scores  associated 
with  the  maximum  independent  VISCAT  I  threat  score  achie 
at  the  fifth  predictor  level 
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Figure  25.  Contingency  table  results  for  the  area  2, 
TAU-48  equal  variance  threshold  model 
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AO  (dependent  data)  and  the  adjusted  VISCAT  I  threat  score 
(independent  data)  for  area  3W,  TAU-24  (PR  model,  MAXPROB  II 
strategy) .  For  this  case  a  grouping  size  of  six  was  selected 
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Figure  28a.  Skill  diagram  and  contingency  table  results  for  area  3W , 

TAU-24  (PR  model)  for  the  MAXPROB  I  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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Figure  28b.  Skill  diagram  and  contingency  table  results  for  area  3W,  TAU-2 
(PR  model)  for  the  MAXPROB  II  strategy.  The  contingency  table 
corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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Figure  28c.  Skill  diagram  and  contingency  table  results  for  area  3W,  TAU-24 
(PR  model)  for  the  natural  regression  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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SELECTION  OP  GROUPING  SIZE 


AO  (dependent  data)  and  the  adjusted  VISCAT  I  threat  score 
(independent  data)  for  area  3W,  TAU-24  (PR+BMD  model,  MAXPROB  II 
strategy) .  For  this  case  a  grouping  size  of  five  was  selected 
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Figure  31b.  Skill  diagram  and  contingency  table  results  for  area  3W,  TAU-24 
(PR+BMD  model)  for  the  MAXPROB  II  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  first 
predictor  level 


contingency  table  corresponds  to  those  scores  associated  with  the 
maximum  independent  VISCAT  I  threat  score  achieved  at  the 
fifth  predictor  level 
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Figure  33.  Contingency  table  results  for  the  area 
TAU-24  equal  variance  threshold  model 
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Figure  34.  Contingency  table  results  for  the  area  3W, 
TAU-24  quadratic  threshold  model 
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Figure  36b.  Skill  diagram  and  contingency  table  results  for  area  3W,  TAU-4 
(PR  model)  for  the  MAXPROB  II  strategy.  The  contingency  table 
corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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Figure  36c.  Skill  diagram  and  contingency  table  results  for  area  3W, 
(PR  model)  for  the  natural  regression  strategy.  The  con 
table  corresponds  to  those  scores  associated  with  the  ma 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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Figure  39a.  Skill  diagram  and  contingency  table  results  for  area  3W,  TAU-48 

(PR+BMD  model)  for  the  MAXPROB  I  strategy.  The  contingency  table 
corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  first 
predictor  level 
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Figure  39b.  Skill  diagram  and  contingency  table  results  for  area  3W,  TAU-48 
(PR+BMD  model)  for  the  MAXPROB  II  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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Figure  39c.  Skill  diagram  and  contingency  table  results  for  area  3W,  TAU-48 
( PR+BMD  model)  for  the  natural  regression  strategy.  The  con¬ 
tingency  table  corresponds  to  those  scores  associated  with  the 
maximum  independent  VISCAT  I  threat  score  achieved  at  the 
fifth  predictor  level 
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Figure  41.  Contingency  table  results  for  the  area  3W, 
TAU-48  equal  variance  threshold  model 
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Figure  42.  Contingency  table  results  for  the  area  3W 
TAU-48  quadratic  threshold  model 
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gure  44a.  Skill  diagram  and  contingency  table  results  for  area  4,  TAU-00 
(PR  model)  for  the  MAXPROB  1  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the 
maximum  independent  VISCAT  I  threat  score  achieved  at  the 
fifth  predictor  level 
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Figure  44b.  Skill  diagram  and  contingency  table  results  for  are  4,  TAU-00 
(PR  model)  for  the  MAXPROB  II  strategy.  The  contingency  table 
corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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Figure  44c.  Skill  diagram  and  contingency  table  results  for  area  4,  TAU-00 

(PR  model)  for  the  natural  regression  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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gure  47a.  Skill  diagram  and  contingency  table  results  for  area  4,  TAU-00 
(PR+BMD  model)  for  the  MAXPROB  I  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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Figure  47b.  Skill  diagram  and  contingency  table  results  for  area  4,  TAU-00 
( PR+BMD  model)  for  the  MAXPROB  II  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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Figure  47c.  Skill  diagram  and  contingency  table  results  for  area  4,  TAU-00 

( PR+BMD  model)  for  the  natural  regression  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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(PR  model)  for  the  MAXPROB  I  strategy.  The  contingency  table 
corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  third 
predictor  level 
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Figure  50b.  Skill  diagram  and  contingency  table  results  for  area  4,  TAU-24 
(PR  model)  for  the  MAXPROB  II  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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Figure  50c.  Skill  diagram  and  contingency  table  results  for  area  4,  TAU-24 

(PR  model)  for  the  natural  regression  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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Figure  51.  Functional  dependence,  AO/Al  statistics  and  96%/05% 
confidence  interval  values  for  area  4,  TAU-24  (PR 
model,  MAXPROB  II  strategy) 
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Figure  53a.  Skill  diagram  and  contingency  table  results  for  area  4,  TAU 
(PR+BMD  model)  for  the  MAXPROB  I  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the 
maximum  independent  VISCAT  I  threat  score  achieved  at  the 
fifth  predictor  level 
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igure  53b.  Skill  diagram  and  contingency  table  results  for  area  4,  TAU- 
(PR+BMD  model)  for  the  MAXPROB  II  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the 
maximum  independent  VISCAT  I  threat  score  achieved  at  the 
fifth  predictor  level 
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Figure  53c.  Skill  diagram  and  contingency  table  results  for  area  4,  TAU 
( PR+BMD  model)  for  the  natural  regression  strategy.  The 
contingency  table  corresponds  to  those  scores  associated 
with  the  maximum  independent  VISCAT  I  threat  score  achieved 
at  the  fifth  predictor  level 
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Figure  55.  Contingency  table  results  for  the  area  4 
TAU-24  equal  variance  threshold  model 
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MRXPROB  II  SCHEME  -  TAU48  -  RRCfl  4 


AO  (dependent  data)  and  the  adjusted  VISCAT  I  threat  score 
(independent  data)  for  area  4,  TAU-48  (PR  model,  MAXPROB  II 
strategy) .  For  this  case  a  grouping  size  of  seven  was  selected 
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Figure  57a.  Skill  diagram  and  contingency  table  results  for  area  4,  TAU-4 
(PR  model)  for  the  MAXPROB  I  strategy.  The  contingency  table 
corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 
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Figure  57c.  Skill  diagram  and  contingency  table  results  for  area  4,  TAU-48 

(PR  model)  for  the  natural  regression  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the  maximum 
independent  VISCAT  I  threat  score  achieved  at  the  fifth 
predictor  level 


Predictor  FD(96)  FD  AO  A0<96)  A1  A1(05) 


NUMBER  OP  PREDICTORS  VS.  SCORING  TECHNIQUES  independent  data 


Skill  diagram  and  contingency  table  results  for  area  4,  TAU 
(PR+BMD  model)  for  the  MAXPROB  I  strategy.  The  contingency 
table  corresponds  to  those  scores  associated  with  the 
maximum  independent  VISCAT  I  threat  score  achieved  at  the 
fifth  predictor  level 
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Figure  60b.  Skill  diagram  and  contingency  table  results  for  area  4,  T, 
( PR+BMD  model)  for  the  MAXPROB  II  strategy.  The  continge 
table  corresponds  to  those  scores  associated  with  the 
maximum  independent  VISCAT  I  threat  score  achieved  at  the 
third  predictor  level 


?  9 


m  r 


10  ^ 
03  00 

6  6 


O  in  vl  m 

<  i- 

<  <  <  < 

SP  jja 

CO  5  K  CO 

9  $  £  -  - 

§  «  d  ?  «? 

«  t—  *  *  n* 

®  ■  ^ 

O  —  m  m  in 


</> 

►- 

in 

m 

00 

< 

< 

< 

T  X 

1  X 

* 

(0 

in 

co 

CO 

00 

O) 

00 

TAU 

wi 

9 

6 

6 

o 

QJ  XJ 

T- 

H 

a 

n 

-X  QJ 

M 

5 

TS1 

10 

h* 

<n 

»- 

X  fn  X  tJ 
(0  QJ 

1167 

104 

15 

n 

Q 

CO 

54 

oo 

UJ 

> 

<N  * 
Ui 

in 

O) 

in 

03 

O 

oi 


iS*D3HOj 


i S*  0  3d  O  i 


(0  H  > 
0)  •  U  OJ 

l-l  >i  O  -H 

<o  O'  w  x 

CL)  C/3  O 


a 

p 

X 

10 

fO 

UJ 

0 

10 

> 

X 

p 

Ul 

QJ 

(X 

X 

0) 

P 

m 

Ul 

U) 

p 

O 

X 

o 

o 

o 

X 

c 

u 

w 

3 

0 

01 

Q)  W  QJ  10 
Li  U)  U)  QJ 
QJ  0  Li 
OJ  Li  X  X 
L  O'*)  |J 
X)  0) 

O  Q  OH 
X  X 
i-i  E-* 

>.  10  Ul  < 

O  liOUH 
C  3  C  CO  QJ 
QJL  OH  > 
Cn  (0  a>  OJ 
C  C  01  X 
•H  Q)  X 
-P  QJ  Li  C  Li 
C  X  Li  QJ  O 
O  X  O  X)  X 
U  U  C  U 
Li  0)  -h 
0  0  4)  CLX> 
C  M-i  i— l  OJ  QJ 
<0  X  'CO  L| 
-«J  C  Q 
e  h  +j  h 
<0  QJ  X 
Li  XJ  >•  6  x 
O'O  U  3LI 

«  E  C  E  H 

■H  OJ  -H  IH 

xj  a  cn  x 

S  C  10  OJ 
-H  CQ  -H  EX 
X  +  X  X 
•H  OS  c  QJ 
JC  A  OX  p 
CO  — -  O  X  10 


sro  wo  s£-q  si  'a  soo-  sz-o- 

l  LUOS  [A  3B03S  IB3HRL  ?  oy  Q3i9nrGB 


Predictor  FD(96)  FD  AO  A0(96)  A1  A1(05) 


interval  values  for  area  4,  TAU-48  (PR+BMD  model, 
MAXPROB  II  strategy) 
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Figure  62 


Contingency  table  results  for  the  area  4 
TAU-48  equal  variance  threshold  model 
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Figure  63.  Contingency  table  results  for  the  area  4, 
TAU-48  quadratic  threshold  model 
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NUMBER  OP  PREDICTORS  VS.  SCORING  TECHNIQUES  independent  data 

TAU2A.1U34  -  AREA  2  -  AAXPROBII-  GROUPING  SIZE  Or  •»  _ 


2  n 

O  T- 


^  to 

Z  *  S  cv 

5  n.  6 

too. 


*“  0) 
ID  CD 


—  M  — 

en  in  i/i 


< 

< 

< 

< 

* 

o 

* 

co 

a> 

CO 

<o% 

io 

CO 

K 

hi 

oo 

o 

K 

o 

it 

o 

o 

H 

<N 

0 

1 

U) 

5 

cn 

CO 

Oi 

h- 

01 

co 

(0 

99 

0) 

CO 

m 

CO 

CO 

0) 

o> 

o 

Ol 

CO 

8 

141 

30 

22 

O 

287 

jSIOIUOj 


isnanoj 


T  ? 


■  * 

v  ; 

•  s,  ^  * 

\  '  : 

\H 


•  O  X 
a  az  x  t— 
2  0  1 
U  <c  H  •  • 

eg  a.  a. 

lj  ■  ■  u  u 

j  1.  CL  O  □ 

U  U  7  Z 

□  a  *—  — • 

iiii 

□  ♦  t  « 


8  0  9'0  >  U  Z  0  G'O  ?'0- 

i  Idas  t  a  3hcjjc,  Ly.jdMi  ?  ay  a3ibnrau 


Figure  64.  Skill  diagram  and  contingency  table  results  for  area  2,  TAU-2 
1984  data  ( PR+BMD  model)  for  the  MAXPROB  II  strategy.  The 
contingency  table  corresponds  to  those  scores  associated  with 
the  maximum  independent  VISCAT  I  threat  score  achieved  at 
the  fourth  predictor  level 
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Figure  68. 


Contingency  table  results  for  the  minimum 
probable  error  threshold  model  (EVAR)  for 
area  4,  TAU-24.  The  contingency  tables  reflect 
VISCAT  I  +  II  versus  VISCAT  III  and  VISCAT  I 
versus  VISCAT  II  separations 
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Figure  69.  Contingency  table  results  for  the  maximum- 
likelihood-  of ■ detection  threshold  model  for 
area  4,  TAU-24.  The  contingency  tables 
reflect  VISCAT  I+II  versus  VISCAT  III  and 
VISCAT  I  versus  VISCAT  II  separations 
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